EPJ manuscript No. 

(will be inserted by the editor) 



00 

g ' Dynamical stability of systems with long-range interactions: 
cn| application of the Nyquist method to the HMF model 

D ' P.H. Chavanis and L. Delfini 

^ I ^ Universite de Toulouse, UPS, Laboratoire de Physique Theorique (IRSAMC), F-31062 Toulouse, France 
'. ^ CNRS, Laboratoire de Physique Theorique (IRSAMC), F-31062 Toulouse, France 
e-mail: chavanis@irsamc.ups-tlse . f r ; Luca.Delf iniOirsamc .ups-tlse . f r 



■ To be included later 

o , 

Abstract. We apply the Nyquist method to the Hamiltonian Mean Field (HMF) model in order to settle 
the linear dynamical stability of a spatially homogeneous distribution function with respect to the Vlasov 
equation. We consider the case of Maxwell (isothermal) and Tsallis (polytropic) distributions and show 
. that the system is stable above a critical kinetic temperature Tc and unstable below it. Then, we con- 

' sider a symmetric double-humped distribution, made of the superposition of two decentered Maxwellians, 

and show the existence of a re-entrant phase in the stability diagram. When we consider an asymmetric 

■ double-humped distribution, the re-entrant phase disappears above a critical value of the asymmetry factor 
' A > 1.09. We also consider the HMF model with a repulsive interaction. In that case, single-humped dis- 

I , tributions are always stable. For asymmetric double-humped distributions, there is a re-entrant phase for 

. 1 < Zi < 25.6, a double re-entrant phase for 25.6 < A < 43.9 and no re-entrant phase for A > 43.9. Finally, 

^ ' we extend our results to arbitrary potentials of interaction and mention the connexion between the HMF 

O , model, Coulombian plasmas and gravitational systems. We discuss the relation between linear dynamical 

O . stability and formal nonlinear dynamical stability and show their equivalence for spatially homogeneous 

distributions. We also provide a criterion of dynamical stability for inhomogeneous systems. 

^ ■ PACS. 5.20.-y Classical statistical mechanics - 05. 45. -a Nonlinear dynamics and chaos - 05.20.Dd Kinetic 

^>f-^ ' theory - 64.60.De Statistical mechanics of model systems 

(N 

CO I 1 Introduction re-entrant phases |55l58j . dynamical phase transitions |60 ) . 

^vq . curious effects due to spatial inhomogeneity |61| . non- 

T— I ■ Systems with long-range interactions have recently been ergodic behaviors |62|31I63(64] , glassy dynamics [6S] , non- 

OO ! the object of considerable interest [112] . These systems ideal effects due to the influence of a thermal bath [29,66- 

■ have applications in different areas of physics, astro- 68] etc... Certainly, the physical richness of these systems 
physics, hydrodynamics and biology Q. Furthermore, they (Pandora box) will lead to further investigations and ap- 
display very peculiar phenomena with respect to ordi- plications. 

^ nary systems with short-range interactions such as neg- 

H : ative specific heats [6-9], inequivalence of statistical en- ^ ^y^^^n^s ^^^^ long-range mteractions, the mean 

. sembles [10-16], numerous types of phase transitions per- ^^^1^ approximation is exact m a proper thermodynamic 

sisting at the thermodynamic limit mM, slow coUisional l^^^* ^,^ere the number of particles N ^ +oo while the 

relaxation [19-25], kinetic blocking due to the absence of strength of the interaction k ~ 1/iV (coupling c onstant ) 

resonances [26-30,25], non-trivial dependence of the col- ^ocs to zero and the volume 1/ ^ 1 remains fixed [112i69j- 

lisional relaxation time with the number N of particles f""' ^ fixed interval of times^ [0, T] and iV ^ +oo, the evo- 

[31-33], long-lived metastable states whose lifetimes in- ^^^lon of the distribution function /(r,v,t) is governed 

crease exponentially with N [SEl], algebraic decay of the (mean field) Vlasov equation. The Vlasov equa- 

the correlation functions |36 l 37 l 28 l 30 l 38) , front structure describes the colhsionless dynamics of the system for 

and slow relaxation of the velocity tails [39130], long-lived sufficiently short times with respect to the colhsional 

quasi stationary states (QSS) [20-23,31,32,40,41,24,25,42], relaxation time: t « tR{N) In fact, the relaxation time 

violent colhsionless relaxation of the Vlasov equation [43- *«(^) mereases rapidly with the nuniber of particles N 

47,29,48-51], out-of-equilibrium phase transitions [52-59], ^^^t' PracticMhe domam of validity of the Vlasov 

equation is huge [69J. Since the Vlasov equation admits 

^ For example, there exists beautiful analogies between self- an infinite number of steady states, the system can be 

gravitating systems, two-dimensional vortices and bacterial found in quasi stationary states (QSS) that differ from 

populations investigated by one of the authors [3-5]. the Boltzmann distribution. These are long-lived stable 
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steady states of the Vlasov equation. For example, if we 
place initially the system in a stable steady state of the 
Vlasov equation /o(r,v), it will remain in that state for 
a very long time until "collisions" (more precisely corre- 
lations due to finite N effects) make it slowly evolve [5^ . 
On the other hand, if we start from an unstable steady 
state of the Vlasov equation, or from an unsteady initial 
condition, the system will spontaneously evolve through 
phase mixing and violent relaxation towards a QSS on 
the coarse-grained scale (see, e.g., |21|3| ). This distribu- 
tion function /qss(''i'^) is a stable steady state of the 
Vlasov equation that can be predicted by the statistical 
theory of Lynden-Bell [l^ if the evolution is ergodic (i.e. 
if the system "mixes" efficiently [69]). As a result, it is 
important on general grounds to study the dynamical sta- 
bility of stationary solutions of the Vlasov equation and 
devise general stability criteria. 

This problem has been first considered in plasma 
physics [70-72] and astrophysics [20-23]. In plasma physics, 
the interaction between like-sign charges is repulsive. How- 
ever, in neutral plasmas, the presence of opposite charges 
screens the interaction between like-sign charges (Debye 
shielding) so that the system is spatially homogeneous. 
Several methods have been developed to study the dy- 
namical stability of spatially homogeneous distribution 
functions /(v) with respect to the Vlasov equation. One 
approach is to determine criteria of nonlinear dynamical 
stability or formal nonlinear stability by minimizing an 
energy-Casimir functional |73j . Another approach is to de- 
termine criteria of linear stability by linearizing the Vlasov 
equation around a steady state, taking the Laplace-Fourier 
transform to obtain the dispersion relation and studying 
the sign of the imaginary part of the complex pulsation 
[70-72] . In some cases, the dispersion relation can be solved 
analytically and the pulsation can be obtained explicitly. 
In more complicated cases, the dispersion relation must 
be solved numerically. We can also have recourse to other 
methods to settle the linear stability of the system with- 
out being required to solve the dispersion relation. In this 
respect, Nyquist [74] has introduced a powerful method to 
investigate the linear stability of a spatially homogeneous 
distribution of the Vlasov equation. This is based on a 
graphical construction that can be easily implemented in 
practice. This method tells whether the distribution is lin- 
early stable (or unstable) but it does not give the value of 
the decay rate (or growth rate) of the perturbation. Using 
this criterion, it can be shown with almost no calcula- 
tion that the Maxwellian distribution is always dynami- 
cally stable in a plasma for perturbations with arbitrary 
wavenumbers. 

For gravitational systems, the interaction between 
masses is attractive leading to possible collapse (Jeans in- 
stability). We could imagine using the Nyquist method 
to settle the dynamical stability of a galaxy described by 
the Vlasov equation. However, we are rapidly confronted 
to the difficulty that a realistic stellar system is spatially 
inhomogeneous [20-23] . This considerably complicates the 
stability analysis and precludes the direct application of 
the Nyquist method. In that case, we must use other meth- 



ods like the Antonov [75-77] criterion of linear dynamical 
stability or methods of nonlinear dynamical stability |78j . 
Using these methods, it can be shown that any distribu- 
tion function of the form / = /(e) with /'(e) < where 
£ = v"^ /2 + ^(r) is the individual energy of a star is lin- 
early [21], and even nonlinearly [78], dynamically stable 
with respect to the Vlasov equation. The Nyquist method 
can be used (see Appendix |D]) to investigate the linear 
dynamical stability of a rather unrealistic infinite and ho- 
mogeneous distribution of stars making the Jeans swindle 
For a Maxwellian distribution, it leads to a graphi- 
cal illustration of the Jeans instability showing that the 
system is unstable with respect to perturbations of suffi- 
ciently low wavenumbers k < kj, where kj = {4'KGf3nY^'^ 
is the Jeans wavenumber. 

In this paper, we consider the Hamiltonian Mean Field 
(HMF) model which is a very much studied toy model 
of systems with long-range interactions |79l29j . In this 
model, the interaction between particles is attractive like 
gravity but, unlike gravity, the system can be spatially 
homogeneous. Therefore, for the HMF model, it is possi- 
ble to use the Nyquist method of plasma physics to in- 
vestigate the linear stability of a homogeneous distribu- 
tion. However, since the interaction is attractive, a crucial 
sign changes in the dispersion relation and this leads to 
new results with respect to plasma physics. As a direct 
application of this method, we can recover by a graphi- 
cal construction the fact that the Maxwellian distribution 
is linearly stable with respect to the Vlasov equation if 
T > Tc and linearly unstable if T < where = ^ 
is a critical temperature. This instability is similar to the 
Jeans instability in astrophysics where the temperature T 
plays the role of the wavenumber k. 

The paper is organized as follows. In Sec. [H we recall 
basic results concerning the Vlasov equation, the disper- 
sion relation and the Nyquist method in relation to the 
HMF model. We first consider the HMF model with an 
attractive interaction (ferromagnetic). Using the Nyquist 
method, we establish general linear stability criteria for 
single-humped and double-humped distributions. In Sees. 
[3] and m we apply the Nyquist method to the Maxwell 
(isothermal) and the Tsallis (polytropic) distributions. We 
show that the system is linearly stable ilT > and lin- 
early unstable if T < Tc where Tc = -tt^ is a critical 
kinetic temperature depending on the polytropic index n 
(the Maxwellian distribution is recovered for n — > -f-oo). 
In Sec. [HI we consider a symmetric double-humped dis- 
tribution made of the superposition of two decentered 
Maxwellians (with separation Va)- We show that this dis- 
tribution has a rich stability diagram in the {va,T) plane 
exhibiting a re-entrant phase. In Sec. [S] we generalize our 
study to the case of an asymmetric double-humped dis- 
tribution where the amplitude of one maximum is larger 
than the other by a factor A > 1. We determine how 
the stability diagram in the (va , T) plane changes as a 
function of A and show that the re-entrant phase disap- 
pears above a very small asymmetry A > Ac = 1.09. In 
Sec. [71 we consider the HMF model with a repulsive inter- 
action (antifcrromagnetic). We show that single- humped 
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distributions (like the Maxwell and the Tsallis distribu- 
tions) are always linearly stable. We also establish general 
linear stability criteria for double-humped distributions. 
More specifically, we investigate the linear stability of a 
double-humped distribution made of two Maxwellian dis- 
tributions. For A = 1 (symmetric case), we show the ex- 
istence of a re-entrant phase. A double re-entrant phase 
appears for A > Zic^-* = 25.6 and they both disappear 
for A > = 43.9. In Sec. E we discuss the exten- 

sion of our results to arbitrary attractive or repulsive po- 
tentials of interaction and make the connexion between 
the HMF model, Coulombian plasmas and gravitational 
systems. We also discuss the relation between linear and 
formal nonlinear stability and show their equivalence for 
spatially homogeneous systems. 



2 Dynamical stability of the HMF model 
2.1 The Vlasov equation 

The HMF model consists in particles of unit mass 
m = 1 moving on a circle and interacting via a cosine 
potential [79129] . The microscopic dynamics of this sys- 
tem is governed by the Hamiltonian equations 



dt 



dH 



dvj 
dt 



dH 



N 9 
V 



cos — 



(1) 



(2) 



where 6i is the angle that particle i makes with an axis 
of reference, Vi = dOi/dt is its velocity and fc > is the 
coupling constant. This system conserves the energy E = 
H and the total mass M = Nm = N. 

The evolution of the A^-body distribution 
Pn{Oi,vi, ...,9N,VN,t) is governed by the Liouville 
equation 



where 

F, = 



dPN 
dt 



N 

E 

i=l 



k 

2^ 



N 



89,. 



F, 



dPi 



N 



dv. 



^sin( 



(3) 



^) (4) 



is the total force experienced by particle i due to the 
interaction with the other particles. From the Liou- 
ville equation, we can construct the BBGKY hierar- 
chy for the reduced distributions Pj(0i, ui, Uj, t) = 

/ PN{di,vi, ...,9N,VN,t)Y\^^j^id9kdvk [2l]- 

We consider the thermodynamic limit A^ — > +oo in 
such a way that the dimensionless energy e = SirE/kAl'^ 
and the dimensionless temperature 77 = kM/AnT are of 
order unity (these scalings are obtained by comparing 



the kinetic and the potential energies in Eq. ([2]) yield- 
ing E ~ N{v'^) ~ NT ~ fcA^^). We can renormahze the 
parameters so that the coupling constant k ~ 1/A^ while 
/3 - 1, i;/A^ - 1 and F = 27r - 1. The dynamical time 
to ~ 27r/ (w^) ~ 1/ ^/kp is also of order unity 0. In the 
thermodynamic limit A^ 4-oo, it can be shown from 
scaling arguments [H] that the A^-body distribution fac- 
torizes in a product of N one-body distributions 



N 



Pn{9i,vi, ...,9N,VN,t) = Y\_Pi{(^i,Vi,t). 



(5) 



i=i 



Substituting this result in the first equation of the 
BBGKY hierarchy, we find that the one-body distribution 
Pi{9,v,t), or equivalently the smooth distribution func- 
tion f{9, V, t) = NPi{9, V, t), satisfies the Vlasov equation 







dt ^' 89 89 8v 



where 



1>{9.t)^- — 



2k: 



<s{9 - 9')p{9',t)d9', 



(6) 



(7) 



is the mean field potential and p{9,t) = / f{9,v,t)dv is 
the spatial density. The mean field average force experi- 
enced by a particle located in 9 is (F) = —^'{9). 

The Vlasov equation admits an infinite number of sta- 
tionary solutions. A spatially homogeneous system {<P ~ 
0) with distribution function (DF) / = f{v) is always a 
steady state of the Vlasov equation. If the system is spa- 
tially inhomogeneous (<? ^ 0), the general form of steady 
distributions of the Vlasov equation is / = /(e) where 
e = v^/2 + <P{9) is the individual energy. In this paper, we 
shall only consider spatially homogeneous distributions so 
that / = f{v). 



2.2 The dispersion relation 

We want to study the linear dynamical stability of a 
spatially homogeneous stationary solution of the Vlasov 
equation described by a distribution function / = /(w). 
Linearizing the Vlasov equation around the steady state 
and taking the Laplace- Fourier transform (writing the per- 



turbation as Sfnuj ~ S(l>nw ^ 
dispersion relation j81|82|29] : 

k 

ein,Lj) = 1 + -((5„,i - <5„ _i) 



'i"S-<^t)^^ we obtain the 



fiv) 



nv 



dv = 0, (8) 



^ In this paper, we use the notations of that are similar to 
those introduced in astrophysics. The link with the notations of 
[32] is made by taking k = 2n/N, 7? = /3/2 and e = 4((7- 1/2). 

^ A rigorous derivation of the Vlasov equation is given by 
Braun & Hepp [50] . 

In fact, we consider here the spectral stability of a distribu- 
tion function / — f{v). For infinite dimensional systems, spec- 
tral stability does not necessarily imply linear stability [73] . 
However, by an abuse of language, we shall often say "linearly 
stable" instead of "spectrally stable". 
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where e{n,uj) is the dielectric function and the integral 
must be performed along the Landau contour (C) |72j . 
There exists solutions only for n = ±1 so that only these 
modes can propagate 0. For the mode n = +1 (n = — 1 
gives the same result), Eq. ([5]) reduces to 



k f f 



dv = 0, 



(9) 



which is the fundamental equation studied in this paper. 
For a given distribution /(w), this equation determines 
the complex pulsation(s) uj = ujr + ii^i of the linearized 
perturbation 5f . Since 5f ~ e"'*, the system is linearly 
stable if loi < and linearly unstable if > 0. For future 
reference, let us recall that 

fm,,f*-m,„, ,„.>„,, (10, 

JcV-UJ J^oo V-UJ 



(11) 



(12) 

2.3 The condition of marginal stability 

The condition of marginal stability corresponds to 

(13) 

where oji is the imaginary part oi lo. In that case, oj = w,, 
is real and the integral in Eq. ^ can be evaluated along 
the real axis by using the Plemelj formula [72] : 



1 



P- 



1 



in6(u — a), 



u — a u — a 
where P denotes the principal value. We obtain 



k 



±SlLdv + nT-f'{uJr)^{). 



(14) 



(15) 



Identifying the real and the imaginary parts, if follows that 
the marginal mode is determined by the equations |29j : 



1 



I'iv) 

V ~ UJr 



f'{uJr) = 0. 



dv = 0, 



(16) 



(17) 



The second relation fixes the pulsation ojr of the perturba- 
tion and the first relation determines the point of marginal 



^ This differs from the case of gaseous systems described by 
the dispersion relation (|160p . 



stability in the series of equilibria 0. Note that the distri- 
bution f{v) can be relatively arbitrary. However, if f(v) 
has a single maximum at w = 0, then ojr = Q (implying 
cij = 0) and the condition of marginal stability becomes 



+ 00 



dv = 0. 



(18) 



2.4 The Nyquist Method 

To determine whether the distribution / = f{v) is stable 
or unstable, one possibility is to solve the dispersion re- 
lation ([9|) and determine the sign of the imaginary part 
of the complex pulsation. This can be done analytically 
in some simple cases [29]. In this paper, we use another 
strategy. We shall apply the Nyquist method introduced 
in plasma physics. This is a graphical method that does 
not require to solve the dispersion relation. The details of 
the method are explained in |72| and we just recall how 
it works in practice. In the e-plane, we plot the Nyquist 
curve [3 (er('^r)j ei(^r)) parameterized by going from 
— oo to -|-oo. This curve is closed and always rotates in 
the counterclockwise sense. If the Nyquist curve does not 
encircle the origin, the system is stable. If the Nyquist 
curve encircles the origin one or more times, the system is 
unstable. The number N of tours around the origin gives 
the number of zeros of in the upper half plane, i.e. 
the number of unstable modes with u)i > Q. The Nyquist 
method by itself does not give the growth rate of the in- 
stability. 



2.5 General properties 

Let us determine some general properties of the Nyquist 
curve for the HMF model. Taking cj^ = 0, we have 

e{iUr) = l + T:P -^-^dv + tn-fiujr). (19) 

2 V-Ur 2 

Therefore, the real and imaginary parts of the dielectric 
function e{ujr) = er('^r) + iu(y>r) are 



k f 



+00 f/ 



fjv) 

V — LOr 



dv. 



ei{uJr) = TT-f'{uJr)- 



(20) 



(21) 



To apply the Nyquist method, we have to plot the curve 
(er(ti^r)j ei(<^r)) parameterized by Ur going from — oo to 
+00. Let us consider the asymptotic behavior for ujr 
±oo. Since f{v) tends to zero for v ±oo, we conclude 
that ei{LOr) ~> for lu^ ±oo and that ei{u!r) > for 



^ When condition (|17|) is fulfilled the integral is convergent 
in V = LJr so that the principal part is not necessary. 
^ This curve is also called an hodograph. 
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ujr — > — oo while ei(wr) < for Ur ^ 
hand, integrating by parts in Eq. (|20 



+00. On the other 
we obtain 



k 



fiv) 



{v — LUr)' 



;dv, 



(22) 



provided that f{v) decreases sufficiently rapidly. There- 
fore, for LOr — > ±00, we obtain at leading order 



r{uJr) 



1 



k p 



{uJr 



±00) 



(23) 



In particular, tr{'^r) ^ 1 for lj^ —>■ ±00. From these re- 
sults, we conclude that the behavior of the Nyquist curve 
close to the limit point (1, 0) is the one represented in Fig. 
[TJ In addition, according to Eq. (PT|) . the Nyquist curve 
crosses the x-axis at each value of utr corresponding to an 
extremum of f{v). For LOr = v^xt, where v^xt is a velocity 
at which the distribution is extremum {f'{vext) = 0), the 
imaginary part of the dielectric function ei{vext) = and 
the real part of the dielectric function 



^{Vext) = 1 + 



fiv) 

' ~ Vext 



dv. 



(24) 



Subtracting the value f {vext) = in the numerator of the 
integrand, and integrating by parts, we obtain 



fjVext) " f{v) 
{V - VextY 



dv. 



(25) 



If vj^iax denotes the velocity corresponding to the global 
maximum of the distribution, we clearly have 



f-r{vMax) < 1. 



(26) 



Finally, for an attractive interaction, we note that if the 
zero of /'(w) corresponds to a maximum (resp. minimum) 
of / then the hodograph crosses the real axis downward 
(resp. upward). On the other hand, if f{v) is symmetric 
with respect to v = 0, then ei((jJr = 0) = 7r|/'(0) = and 



e',{cjr - 0) = 



+00 f'(v) 



dv = SO that the Nyquist 



curve has a vertical tangent at (6^(0), 0) when ujr — 
provided that e^(cj^ = 0) cx /"(O) ^ 0. 



2.6 Single-humped distributions 

Let us assume that the distribution f{v) has a single max- 
imum at V = Vo (so that f {vq) = 0) and tends to zero for 
V — > ±00. Then, the Nyquist curve cuts the x-axis {eiiuor) 
vanishes) at the limit point (1,0) when lo^ — > ±00 and 
at the point (er(fo)i 0) when ujr = vq. Due to its behavior 
close to the limit point (1,0), the fact that it rotates in the 
counterclockwise sense, and the property that £^(1^0) < 1 
according to Eq. (pS)) . the Nyquist curve must necessarily 
behave like in Fig. [TJ Therefore, the Nyquist curve starts 
on the real axis at er(^r) = 1 for ujr — > —00, then going in 
counterclockwise sense it crosses the real axis at the point 
er(wo) < 1 and returns on the real axis at er(wr) = 1 foi' 



LOr +00. According to the Nyquist criterion exposed in 
Sec. 12.41 wc conclude that a singlc-humpcd distribution is 
linearly stable if 



er(i'o) 



+00 



V - Vq 



dv > 0, 



(27) 



and linearly unstable if £^(^^0) < 0. The equality corre- 
sponds to the marginal stability condition p6 |) -(fT7 ll . When 
the system is unstable, there is only one mode of instabil- 
ity since the Nyquist curve rotates only one time around 
the origin. This stability criterion, which is valid for single- 
humped distributions, was stated without proof in [29| . It 
is here justified from the Nyquist method. In particular, a 
symmetric distribution / — f(v) with a single maximum 
at Wo = is linearly dynamically stable iff 



+00 f/ 



fiv) 



dv > 0. 



(28) 



This coincides with the criterion of formal stability given 
in |32|29j . Therefore, for spatially homogeneous distribu- 
tions, the criterion of linear dynamical stability coincides 
with the criterion of formal nonlinear dynamical stability 
(see SecE]). 



2.7 Double-humped distributions 

Let us consider a double-humped distribution with a 
global maximum at VMax , a minimum at Vmin and a local 
maximum at Vmax- We assume that VMax < Vmm < Vmax- 
The Nyquist curve will cut the a;-axis at the limit point 
(1,0) and at three other points {^rivMax),^), {^r(vmin),Q) 
and (er(wmaa;), 0). Wc also know that the Nyquist curve 
can only rotate in the counterclockwise sense and that 
^r{vMax) < 1 according to Eq. (|26p . Then, we can convince 
ourselves by making drawings of the following results. If 

(+ + +): ^rivMax) > 0, er{v„iin) > 0, er{Vrnax) > 0, 
(,H )'■ (-r\VMax ) > 0, er{Vmin) < 0, er{v,nax) < 0, 

( h): Cr^VMax) < 0, eriv„iin) < 0, eriVrnax) > 0, 

(H |-)E|: er{vMax) > 0, er{Vmzn) < 0, er{v 

max ) > 0, 

the Nyquist curve does not encircle the origin so the sys- 
tem is stable. If 

( ): ^rivMax) < 0, er{Vmin) < 0, CriVmax) < 0, 

(- + +): ^rivMax) < 0, eriVmin) > 0, eriVmax) > 0, 
(+ + -): <^r{vMax) > 0, er{Vmin) > 0, Cr^V 

max ) < 0, 

the Nyquist curve rotates one time around the origin so 
that there is one mode of instability. Finally, if 

(- + -): er{vMax) < 0, er{Vmin) > 0, er{Vmax) < 0, 

the Nyquist curve rotates two times around the origin so 
that there are two modes of instability. Cases (-1- + +), 
( ); ( — h +) and ( — \ — ) are observed in Sec. [Slfor 



This case is relatively tricky. It is obtained by making two 
loops in the S-E and N-W quadrants of the e-plane (assuming 
that this is possible for some double- humped distributions). 
The origin is then encircled by a sort of triangle but since the 
curve rotates clockwise along this triangle, the origin is in fact 
exterior to the Nyquist curve. Thus the system is stable. 
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an asymmetric double-humped distribution made of two 
Maxwellians. The other cases cannot be obtained from 
this distribution but they may be obtained from other 
distributions. 

If the double-humped distribution is symmetric with 
respect to the origin with two maxima at and a min- 
imum at V = 0, only three cases can arise. If 

(-h-H-H): er{v^) > 0, er{0) > 0, 

(+-+): e,(«,) >0, e,(0) <0, 
the Nyquist curve does not encircle the origin so the sys- 
tem is stable. If 

( ): < 0, 6,(0) < 0, 

the Nyquist curve rotates one time around the origin so 
that there is one mode of instability. Finally, if 

{- + -): erM<0, e,(0)>0, 
the Nyquist curve rotates two times around the origin so 
that there are two modes of instability. Cases (+ + +), 

( ) and ( — I — ) are observed in Scc.[5]for a symmetric 

double-humped distribution made of two Maxwellians. 



2.8 A particular solution of e{uj) = 

For the attractive HMF model we can look for a solution of 
the equation e(cij) = in the form lu = itOi corresponding 
to LOr = 0. In that case, the perturbation grows [ui > 0) 
or decays {uji < 0) without oscillating. For uji > 0, the 
equation e{uj) ~ becomes 



1 + - / ^ ^ dv^O. 
2 V - luji 



(29) 



Multiplying the numerator by v + iuji and separating real 
and imaginary parts, we obtain 



k f+°° vf(v) , 

, + Ujf 



fiv) 



dv = 0. 



(30) 



(31) 



If we consider distribution functions f{v) that are sym- 
metric with respect to w = 0, Eq. pTjl is always satisfied. 
Then, the growth rate uJi is given by Eq. ([50)1 . 
For ijji <Q. the equation e{u)) = becomes 

1 + -/ dv + ikTTf'{iu;,) = 0. (32) 



Multiplying the numerator by v + iuJi, assuming that f{v) 
is even, and separating real and imaginary parts, we obtain 



Note that the solution uj ~ iuji exists only for attrac- 
tive interactions. Note also that there may exist other 
solutions to the equation e(a;) = 0. However, for unsta- 
ble single-humped distributions, the Nyquist curve encir- 
cles the origin only once (see Sec. 12. 6p so that w = iuJi 
with > is the only solution of e{uo) = (for stable 
single-humped distributions there may be other solutions 
of e(a;) = with ^ and LOi < 0). Explicit solutions of 
e(ti') = with UJ = iuji are given in [29] for the Maxwell 
and Tsallis distributions. 



3 The Maxwell distribution 

We consider the Maxwell (or isothermal) distribution 

/ O \ 1/2 

/(«)=Uz P^~^"^^ (34) 



where p = A//(27r) is the spatial density and T = = 
{v^) is the kinetic temperature 0. We justify here the 
Maxwell distribution as a particular stationary solution of 
the Vlasov equation [29|83j . The Maxwell distribution rep- 
resents also the statistical equilibrium state of the system. 
In that case, T can be interpreted as the thermodynamical 
temperature. 

The Maxwellian distribution has a single maximum 
at w = 0. Therefore, the condition of marginal stability 
([T7)) implies LUr = 0. From Eq. PB)) . we find that the 
Maxwellian distribution is marginally stable for T = 
where we have introduced the critical temperature 



kM 
47r ' 



(35) 



According to the criterion ([^5]) . the Maxwell distribution 
is linearly dynamically stable if T > Tc and linearly dy- 
namically unstable if T < Tc. 

The dielectric function ^ associated to the 
Maxwellian distribution is 



, . 1 k ( P 



1/2 



-^e-^^ dv. (36) 



Introducing the dimensionless temperature 

l3kM 



11 



47r ' 



and the dimensionless pulsation and velocity 



[2 



/ 47r 



1/2 



V 



/ 47r 
[kM 



1/2 



(37) 



(38) 



1 + ^ 



+ 00 



vfiv) 
2J_^ v^+tof 



For any distribution function /(v), we define the kinetic 
dv + ikirfiioji) = 0, (33) temperature by the relation im{(v - (v))^) = ^Tki„ where 

d is the dimension of space. The energy (conserved control 
parameter) is then E = ^Nm{v^) = ^Nm{vf + ^NTkin- 



and Eq. The equation for the imaginary part is al- 

ways satisfied. Then, the damping rate tUi is given by Eq. 



For one dimensional systems (d = 1) and particles of unit 
mass (m =1), the kinetic temperature is equal to the velocity 
dispersion Tkin ~ {v^) when {v) — 0. 
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Fig. 1. Nyquist curve for the Maxwellian distribution (|34p . 
The DF is stable for rj < rjc {T > Tc), marginally stable for 
ri = rjc {T = Tc) and unstable for rj > rjc {T < Tc). We have 
taken r; = 2, 1, 0.5 from the outer to the inner curve. 



the dielectric function can be rewritten 

eif2) = 1 - l^wi^^^2), 

where 

1 f X 



W{z) 



111 J c — ^ 



-e 2 dx 



(39) 



(40) 



is the M^- function of plasma physics [70l71j . We note that 
VF(0) = 1. For any complex number z, we have the ana- 
lytical expression 



W{2 



1- ze 



2 dx + i\ —ze 2 



(41) 



When Hi = 0, the real and imaginary parts of the 
dielectric function e{f2j.) = er{f2r)+i^i{^r) can be written 



with 



e^{Qr) = -r]W^{y/rjnr), 



Wr{z) ^ 1 - ze 2 / e 2 dx, 
Jo 



W,iz) 



(42) 
(43) 

(44) 

(45) 



where z is here a real number. The condition of marginal 
stability corresponds to er{f^r) = ti{^r) = 0. The condi- 
tion ei{f2r) = 0, which is equivalent to f'{Qr) ~ 0, implies 
f2r = 0. Then, the condition er(i7,.) = leads to r] = ijc 
with 

Vc = l- (46) 

To apply the Nyquist method, we need to plot the 
curve {€r{f^r),^i{^r)) in the e-plane. For 12^ ±oo, this 
curve tends to the point (1,0) in the manner described in 
Sec. 12.61 On the other hand, for fir — 0, it crosses the x- 
axis at (er(0) = 1 — 77, 0). The Nyquist curve is represented 
in Fig. [1] for several values of the inverse temperature rj. 



For 7] < 1 (i.e. T > Tc), the Nyquist curve docs not encircle 
the origin so that the Maxwellian distribution is stable. 
For ?7 > 1 (i.e. T < Tc) the Nyquist curve encircles the 
origin so that the Maxwellian distribution is unstable. For 
7] — 1 (i.e. T = Tc) the Nyquist curve passes through the 
origin so that the Maxwellian distribution is marginally 
stable. 

The stability limits obtained with the Nyquist method 
for isothermal and polytropic (see next section) distri- 
butions coincide with those found in [29j . However, the 
Nyquist method docs not give the value of the growth 
rate (in the unstable regime) or decay rate (in the sta- 
ble regime) of the perturbation. These values have been 
obtained in [29] by solving the dispersion relation ([9]) . 



4 The Tsallis distributions 

We consider the Tsallis (or polytropic) distributions writ- 
ten in the form (see [IS] and Appendix IB]) : 



fiv) = Br, 



V2ttT 



2(n+ l)r_ 



1-1/2 



(47) 



where P = is the density, T = 1/(3 ~ (v^) is the kinetic 
temperature and _B„ is a normalization constant given by 



B„ 



Bn — 



r{n + i) 



r(n + l/2)(n -1-1)1/2' 

r(l/2- n) 
r(-n)[-(n+l)]i/2' 



1 

n > — , 
2 



n < -1. 



(48) 
(49) 



We limit ourselves to these indices so that the distribu- 
tions arc normalizable and their second moments (w^) ex- 
ist. Some representative distributions are plotted in Figs. 
[5J[3]and|3J For n > 1/2 the distributions have a compact 
support since they vanish at u = ±Vm where 



y/2in+l)T. 



(50) 



For 1) > Vm, we set / = 0. Therefore, the notation [x]+ = x 
for a; > and [x]+ = for a; < 0. For n < —1, the dis- 
tribution functions remain strictly positive for all v and 
they decrease algebraically like jwl"*^!"^") for v — > ±00. 
For n ±00, we recover the Maxwell distribution 
We justify here the Tsallis distributions (|47]l as partic- 
ular stationary solutions of the Vlasov equation |29|83j . 
These distributions possess a lot of interesting mathemat- 
ical p rop erties so that their study is interesting in its own 
right [3- We also note that the Tsallis distributions (jTT)) 



^'^ Such distributions have been extensively studied in astro- 
physics because they represent an important class of stationary 
solutions of the Vlasov- Poisson system called stellar polytropes 
[21j . Stellar polytropes have been introduced by Plummer [84] 
and Eddington 05]. However, these distributions were writ- 
ten in a form that is less convenient than the "Tsallis form" 
(|215|l or (|47|. The nice forms (|215|l or ([47]) allow us to make 
a smooth connexion between polytropic and isothermal distri- 
butions for g ^ 1 or n — > 00. See [83] for further discussion on 
these historical issues. 
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Fig. 2. Tsallis distribution with n < — 1 (power-law tail) and 
Maxwell distribution (corresponding to n = ±oo). 



0.4 



0.3 



iO.2 



0.1 



1 ' 1 ' 1 
n > 3/2 
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^ I.I.I 


1 . 1 V. 



Fig. 3. Tsallis distribution with n > 3/2 (specifically n = 3). It 
has a compact support with a horizontal tangent at v — ±v,n. 
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1 1 1 


■ -V__^ / 
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Fig. 4. Tsallis distribution with 1/2 < n < 3/2 (specifically 
n = 1.2). It has a compact support with a vertical tangent at 

V = ±Vm. 



have been observed in numerical simulations of the HMF 
model either in case of incomplete relaxation [31|33j or 
during the coUisional evolution of the system [60] (see dis- 
cussion in Sec. 3.4 of [24|). 

The Tsallis distributions have a single maximum at 
u = 0. Therefore, the condition of marginal stability (|17p 
implies oj^ = (it is shown below that the minimum at 



V = ±Vm for n > 3/2 docs not lead to marginal stability). 
From Eq. (fT6|) . we find that the Tsallis distribution of 
index n is marginally stable for T = where we have 
introduced the critical temperature 



kM 



1 47r 



(51) 



For n ±cxD, we recover the critical temperature psp 
corresponding to the Maxwell distribution. On the other 
hand, according to the criterion ([SS]) . the Tsallis distribu- 
tion of index n is linearly dynamically stable if T > Tc 
and linearly dynamically unstable if T < Tc. 

The dielectric function associated to the Tsallis distri- 
bution of index n is 



e{uj) = 1 - 



zBr, n 



1 



1 



2(n+l)T 



2 j (n + l)r 

n-3/2 

+ 



■ dv. 



(52) 



Introducing the dimensionless temperature (|37p and the 
dimensionless pulsation ([38]) it can be rewritten 



e(/2) = 1 



where 



1 Bn / 1 

— [ n — - 
2t: n \ 2 



a;[l 



1-3/2 



2(n+l)J + 



X — Z 



(53) 

-dx, 
(54) 



generalizes the W function of plasma physics to the case 
of Tsallis distributions We note that M^(")(0) = 1. 

When i7j = 0, the real and imaginary parts of the 
dielectric function e{f2r) = €r{^r) + i^i{^r) can be written 



with 

W,(")(-z) = 



/27r n 



1 



n - 
n 



n + 



(55) 
(56) 



P 



iri-3/2 



2(ri+l) 



-dx, 
(57) 



(")/-,^ - 



2(n+ 1) 



-3/2 



(58) 

where z is here a real number. The condition of marginal 
stability corresponds to er{f2r) = ^ii^r) = 0. The condi- 
tion ei(f]r) = 0, which is equivalent to f {Qr) = 0, implies 
fl,. = 0. Then, the relation er{fir) = leads to = rjc with 



1 



(59) 
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-1 -0.5 0.5 1 1.5 

£^ 

Fig. 5. Nyquist curve for the Tsallis distributions with n < — 1 
(specifically n = —2 yielding rjc — 1/2). The DF is stable for 
r] < rjc, marginally stable for rj = r]c and unstable for rj > rjc- 
We have taken 77 = 1, 0.5, 0.4 from the outer to the inner curve. 



-0.5 0.5 1 1.5 2 

e 

r 

Fig. 6. Nyquist curve for the Tsallis distribution with n > 5/2 
(specifically n = 5 yielding rjc = 6/5). The DF is stable for 
rj < rjc, marginally stable for rj = rjc and unstable for rj > rjc. 
We have taken 77 = 2, 6/5, 0.5 from the outer to the inner curve. 
The curve has a horizontal tangent at {er{Vm),0)- 



To apply the Nyquist method, we need to plot the 
curve (eri^r), £ii^r)) in the e-plane. We have to distin- 
guish different cases: 

• For n < — 1, the Tsallis distributions have a single 
maximum at ij = and they decrease algebraically at 
infinity. For [2,. ±00, the curve {er{f2r),ii{^r)) tends 
to the point (1,0) in the manner described in Sec. 12.61 
On the other hand, for fi^ = 0, it crosses the cc-axis at 
(er(0) = 1 — rj/rjc,0). The Nyquist curve is represented 
in Fig. [5] for several values of the inverse temperature t]. 
For rj < rjc (i.e. T > Tc), the Nyquist curve does not 
encircle the origin so that the Tsallis distribution is stable. 
For rj > rjc (i.e. T < Tc) the Nyquist curve encircles the 
origin so that the Tsallis distribution is unstable. For 77 = 
Tjc (i.e. T = Tc) the Nyquist curve passes through the 
origin so that the Tsallis distribution is marginally stable. 
Therefore, the Nyquist curve of polytropic distributions 
with index n < —1 is very similar to the Nyquist curve of 
Maxwellian distributions (recovered for n = —00). They 
essentially differ in the value of the critical temperature 
'7c- 

• For 71 > 3/2, the Tsallis distributions have a single 
maximum at v ~ and they vanish at t; = ±Wm with 
a horizontal tangent: f'{ztvm) = (sec Fig. |3]). There- 
fore, when i7,. — *■ ±Vm, the Nyquist curve (e,.(/2r), ^ii^r)) 
tends to the point {er{Vm) = 1 + iv/Vc)Cn,0) where 
Cn = -w!:"\y/2{n + 1)) > 0. This parameter is plot- 
ted as a function of n in Fig. [^5] of Appendix [Ul We note 
that (]?,.) > for Qr and €^(^7^) < for 
^2r — + Vj^. We need now to distinguish three sub-cases. 
For n > 5/2, f"{vm) = 0: this implies e'-{Vm) = so that 
the Nyquist curve has a horizontal tangent at {er{Vm), 0). 
For 3/2 < n< 5/2, /"(vm) = +00: this implies e'^{V„t) = 
00, so that the Nyquist curve has a vertical tangent at 
(e,(F„),0). For n = 5/2, f"{v,n) = 15p/ {Uy/lT^^^) and 
e'.(F„) = 157r7/3/2/(i4^) are finite so that the Nyquist 

We note that er(Kn) > 1 (it is equal to unity for n 00) 
so that fir — ±Vm can never correspond to a marginal mode 
(it satisfies ei = but it can never satisfy tr = 0). 
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Fig. 7. Nyquist curve for the Tsallis distribution with 3/2 < 
71 < 5/2 (specifically n = 2 yielding rjc = 3/2). The DF is 
stable for rj < rjc, marginally stable for rj = rjc and unstable 
for rj > rjc- We have taken rj = 2, 3/2, 0.5 from the outer to the 
inner curve. The curve has a vertical tangent at (er(Kn),0). 



0.1 - 



oj' 0- 



-0.1 - 




1.2 1.4 1.6 1.8 

Fig. 8. Behavior of the Nyquist curve close to the point 
(er(Vm),0). For n > 5/2 (specifically n — 3) the tangent is 
horizontal, for n — 5/2 the tangent has a finite slope and for 
3/2<n<5/2 (specifically n = 2) the tangent is vertical. 
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Fig. 9. Nyquist curve for the Tsallis distribution with n = 3/2 
yielding rjc = 5/3. The DF is stable for rj < rjc, marginally 
stable for r] = rjc and unstable for r] > rjc- We have taken 
r? = 2, 5/3, 0.5 from the outer to the inner curve. 



curve has a finite slope at (6^(1^™), 0). These three behav- 
iors are represented in Figs.[6l[7]and[8l On the other hand, 
for Qr < —Vrn, the imaginary part of the dielectric func- 
tion ei{f2r) = and the real part er{fir) varies from 1 (for 

fir -OO) to er{Vm) = 1 + {v/Vc)Cn (for f2r -V~); 

see Appendix [C] For i?^ > Km , the imaginary part of the 
dielectric function e,;(]7r) = and the real part eri^^r) 
varies from 6^(1/™) = 1 + (r?/?7c)Cn (for J?^ ^ to 
1 (for f2r — > +oo). This corresponds to a segment on 
the X-axis betvifeen (1,0) and (1 -I- ('7/'7c)Cn: 0) iii Figs. 
[5] and [7] (for clarity, we have slightly shifted the segment 
from the x-axis in Fig. [7] to show the complete path). 
Finally, considering the value fi,. = 0, we see that the 
Nyquist curve crosses the x-axis at (er(0) = 1 — rj/rjcO). 
For r] < r]c (i.e. T > Tc), the Nyquist curve does not 
encircle the origin so that the Tsallis distribution is sta- 
ble. For ij > rjc (i.e. T < Tc) the Nyquist curve encircles 
the origin so that the Tsallis distribution is unstable. For 
Tj = rjc (i.e. T ~ Tc) the Nyquist curve passes through 
the origin so that the Tsallis distribution is marginally 
stable. The Nyquist curve of polytropic distributions with 
index n > 5/2 is very similar to the Nyquist curve of 
Maxwellian distributions (recovered for n = -l-oo). They 
essentially differ in the presence of a segment between 
(1,0) and (1 -|- {r]/r]c)Cn,0) and in the value of the crit- 
ical temperature rjc- For n < 5/2, the Nyquist curves of 
polytropic and isothermal distributions become sensibly 
different. 

• For 1/2 < ?7 < 3/2, the Tsallis distributions have a 
single maximum at w = and they vanish at v = 
with a vertical tangent: fi^Vm) = ±oo (see Fig. HJ for 
n = 3/2, the slope /'(T^m) = ±3p/(10r) is finite). 
Therefore, the imaginary part of the dielectric function 
ei{f2r) — ^ ±oo when fi^ — ^ TKn (for n = 3/2, the imag- 
inary part of the dielectric function ei{f2r) ±37r77/10 
when Sir TVvi)- We need now to distinguish three 
subcases. For 1 < n < 3/2, the real part of the di- 
electric function er{Qr) +oo when f2r ±Kn- For 
n ~ 1, the dielectric function er{f2r) is independent on 
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Fig. 10. Nyquist curve for the Tsallis distribution with 1 < 
n < 3/2 (specifically n = 1.2). The DF is stable for rj < rjc, 
marginally stable for rj — rjc and unstable for rj > rjc. 
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Fig. 11. Nyquist curve for the Tsallis distribution with n — 1 
yielding rjc = 2. The DF is stable for rj < rjc, marginally stable 
for rj = rjc and unstable for rj > rjc- We have taken 77 = 3, 2, 0.5 
from the outer to the inner curve. 
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Fig. 12. Nyquist curve for the Tsallis distribution with 1/2 < 
n < 1 (specifically n = 0.7). The DF is stable for rj < rjc, 
marginally stable for rj = rjc and unstable for rj > rjc. 
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Fig. 13. Nyquist curve for the Tsallis distribution with ?i = 
1/2 yielding rjc = 3. The DF is stable for -q < rjc, marginally 
stable for rj = r]c and unstable for rj > rjc- 



G [—Vm,Vm]- For 1/2 < 71 < 1, the dielectric func- 
tion er{f2r) —oo when f2r iVm', see Appendix [Cl 
This leads to the Nyquist curves represented in Figs. [9l 
[ini [ni and [TH On the other hand, for ilr < -Vm, the 
imaginary part of the dielectric function ei{f2r) = and 
the real part eri^r) varies from 1 (for f2r —oo) to 
+00 (for f2r —Vj^). For f2r > Kn, the imaginary 
part of the dielectric function ei(J7r) = and the real 
part er{f2r) varies from +00 (for Hr — > +Vj[) to 1 (for 
i7r — > +00). This corresponds to a segment on the a;-axis 
between (1,0) and (+00, 0) in Figs. El [lOl [U] and [l2] (for 
clarity, we have slightly shifted the segment from the x- 
axis to show the complete path). Finally, considering the 
value Hr = 0, we see that the Nyquist curve crosses the 
X-axis at (er(0) = 1 - ri/r]c,0). For r] < r]c (i.e. T > Tc), 
the Nyquist curve does not encircle the origin so that the 
Tsallis distribution is stable. For 77 > 77c (i-e. T < Tc) the 
Nyquist curve encircles the origin so that the Tsallis dis- 
tribution is unstable. For rj = rjc (i.e. T = Tc) the Nyquist 
curve passes through the origin so that the Tsallis distri- 
bution is marginally stable. Note that it is important to 
determine the complete path of evolution, from f2r — —00 
to f2r = +00, in order to apply the Nyquist criterion. 

• The index n ~ 1/2 deserves a particular attention. 
In that case, / = 770 for l^'l < Vm and / = other- 
wise. It corresponds to the Fermi distribution or to the 
water-bag model. The constant 770 is determined by the 
normalization condition p = J f dv yielding p = 2r]ov„ 



The kinetic temperature T 



is given by T = w„j/3. 



The derivative of the distribution function is f'{v) = 
rjo [d{v + vq) — 6{v — vo)]. According to Eq. (fTS]) . the crit- 
ical temperature is Tc = kM/12Tr i.e. 77^ = 3 in complete 
agreement with Eqs. ([5T|) and ifSQ]) . More precisely, for 
n = 1/2, we can obtain an analytical expression of the 
dielectric function in the form 1291: 



with 



(60) 
(61) 



The condition e{uj) = determines the complex pulsation. 
For T > Tc, the system is stable and 

Lj ±V3iT -Tc)^^^. (62) 
For T < Tc, the system is unstable and 

tu = ±iV3{Tc - T)i/2. (63) 
On the other hand, for f2i = 0, we get 



1 1 



3 ' 1 - hjf2^ 



e^iHr) = 0. (64) 



Therefore, the Nyquist curve is made of two segments 
00, 1 — 77/3] and [1, +oo[ as represented in Fig. 1131 



5 The symmetric double-humped distribution 
5.1 Determination of the extrema 

We consider a symmetric double-humped distribution of 
the form 




(65) 



It corresponds to the superposition of two Maxwellian dis- 
tributions with temperature T — l/f3 centered in Va and 
—Va respectively (see Fig. [M]). This distribution models 
two streams of particles in opposite direction. The av- 
erage velocity is (v) = and the kinetic temperature 
Tkin = {v'^) = T + v^. The velocities vq at which the 
distribution function f{v) is extremum satisfy /'(vo) = 0. 
They are determined by the equation 



Va+Vo' 



(66) 



Introducing the dimensionless temperature (|37[) . the di- 
mensionless velocity ()38|) and the dimensionless separa- 
tion 



a = Va 



/ 47r 
\kM 



1/2 



Eq. (|66p can be rewritten 
1 

^ ^ 2aVn 



In 



a + Vo 
a-Va 



(67) 



(68) 



For a given value of the inverse temperature 77 and sepa- 
ration a, this equation determines the velocities Vq where 
f{V) is extremum. We note that Vq G [—a, +a]. It is con- 
venient to introduce the variables 



x = Vo/a, y^i]a^ 



(69) 



For fixed a, the parameter y plays the role of the inverse 
temperature and the parameter x plays the role of the 
velocity. Then, we have to study the function 



1 (l + x 
2x \1 — X 



(70) 
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Fig. 14. Symmetric double-humped distribution made of two 
Maxwellian with separation a. If r; > the distribution 

has two maxima at ±V, and one minimum at Vo = while for 
77 < it has only one maximum at Vb = 0. 



for a; e] — 1, +1[. This function is plotted in Fig.[T5l It has 
the following properties: 



y{~x) = y{x), 



y{x) iln(l 



r 



(71) 

(72) 
(73) 



y(o) = 1. 

The extrema of the distribution function (|65p can be deter- 
mined from the study of this function. First, considering 
Eq. we note that f{V) always has an extremum at 

Vg = 0, for any value of rj and a. This is a "degenerate" so- 
lution of Eq. ([70)) corresponding to the vertical line a; = 
in Fig. [151 On the other hand, if y > 1, i.e 77 > 1/a^, there 
exists two other extrema Vq = ±1^ where 14 = ax* with 

In conclusion, for a given separation a: 

• ii r] > the distribution function f{V) has two 
maxima at Vq = ±14 and one minimum at Vq = 0. 

• if 77 < 1/a^, the distribution function f{V) has only 
one maximum at Vq = (the limit case rj = Ija? corre- 
sponds to /"(O) = 0). 

Remark: We can note the formal analogy with the 
mean field theory (ferromagnetic transition) of the one 
dimensional Ising model without magnetic field (compare 
Eq. ([Ml) with Eq. (14) of |86]). 



5.2 The condition of marginal stability 

The dielectric function associated to the symmetric 
double-humped distribution ([65p is 

e{f2) = 1 - I [W{^j{Q - a)) + W{^{n + a))] , (74) 

where W{z) is defined in Eq. ([40|) . When Qi = 0, the real 
and imaginary parts of the dielectric function e{Qr) = 
Cri^r) + iei{Qr) Can be written 



er{f2r) = l-j [Wr{^{[2r - a)) + Wr{^{Qr + «))] 
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Fig. 15. The function y{x) for the symmetric double-humped 
distribution. 



e,{f2r) = [W,(V^(r2, - a)) + W,{^{f2r + a))] , 

(76) 

where Wr{z) and Wi{z) are defined in Eqs. ([44|) - ([45|) where 
z is here a real number. The condition of marginal stabil- 
ity corresponds to er{f^r) = ^ii^r) = 0. The condition 
£i{f2r) = is equivalent to 



The condition eri^r) ~ leads to 



(77) 



1 - j [Wri^jif2r. - a)) + Wri^ifir + a))] = 0. (78) 



Therefore, according to Eq. ([77p . the real pulsation Qr 
is equal to a velocity Vq at which the distribution ([S5)) 
is extremum. The second equation (|78p determines the 
value(s) r]c{a) of the temperature at which the distribution 
is marginally stable. 



5.2.1 The case uj^ ~ 

Let us first consider the value [2,- = that is solution of 
Eq. ([77]) for any a and r/. In that case, Eq. ([78]) becomes 



vWriVva) = 1, 



(79) 



where we have used Eq. (|207p . For fixed a, this equation 

determines the temperature(s) r]i^\a) corresponding to 
a mode of marginal stability with fir = 0. Introducing 
y = rfo? > 0, the foregoing equation can be rewritten in 
the parametric form 



a^^yWriVy), 



WriVy)' 



(80) 
(81) 



The function defined by Eq. ([80|) is plotted in Fig. [T6| 



It vanishes at y 



2/r, 



where = 1.307 is the 



(75) 



zero of Wr{z) (see Appendix [X)) . For y > ymax = 1.708, 
the function ([80]) is negative so we restrict ourselves to 
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Fig. 16. — yWr{y/y) as a function of y = rjo? 
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Fig. 17. »7 — l/Wri^/y) ^ ^ function ol y — i]a^ 



y G [0,2/maa:]- In this Interval, it has a maximum at 
y = Vm = 0.724 and its value is qm = 0.556. Therefore, 
for a > Qm, Eq- ((80|) has no solution. For a = ajv/, Eq. (|80p 
has one solution y ~ yM determining through Eq. (|8ip a 
temperature 77^"'**^^ = 2.338. For a < qm, Eq. (l80|) has 
two solutions yi^2- Then, Eq. (j8ip determines two temper- 
atures r]c^^^'^'^\a) as illustrated in Fig. [ITl These are the 
solutions of Eq. (|79p for a < um . 

Let us consider the limit a —^ 0. In that case, we see 
in Fig. 1161 that yi — s- and 2/2 ~* ymax- In the first case 
j/i 0, using Eq. ((2T4)) . we find from Eq. ([81]) that 



^7, 



(o)(i) 



(a) 



1, 



(a^O). 



(82) 



This result is to be expected since, for a = 0, the distribu- 
tion ([55)) reduces to the Maxwellian. We thus recover the 
critical temperature (|46p . In the second case 1/2 — > ymax, 
we find from Eq. ([80]) that ~ ymaxWri^/y-i) ■ Sub- 
stituting this equivalent in Eq. (|8ip . and recalling that 



we get 



J0)(2) 



(a^O). 



(83) 



For a 0, this solution rjc"''"'' (a) — > -fcx) so that it does 
not appear in the study of the Maxwellian distribution 
corresponding to a = 0. 



In conclusion: 

• if a < om, there exists two critical temperatures 
T]i°^'''^\a) and 77^°^^^^ (a) determined by Eqs. (l80l)-(l8T]) cor- 
responding to a marginal mode (i?,, = 0, i?^ = 0). These 

two branches merge at the point (aj\/, r^c"''^*^^). 

• if a > aj\f , there is no marginal mode {f2r — 0, fit ~ 

0). 



The curves jyi'^''^^'' (a) and •q^'^'^'ia) corresponding to 
the marginal mode with zero pulsation /2r = are plotted 
in Fig. [HI 



(0)(2) 



5.2.2 The case cj^ 7^ 

We now consider the pulsations Qr ~ ±K that are solu- 
tions of Eq. (j77|) for 77 > To determine the tempera- 
ture(s) at which the distribution (|65|) is marginally stable, 
we have to solve 

1 - \ [Wr{V^[V, - a)) + Wr[^{V, + a))] = 0, (84) 



where K is given by 



1 



2aK 



■In 



(85) 



Eliminating between these two expressions, we obtain 
the critical temperature (s) 'q^\a) as a function of a. How- 
ever, it is easier to proceed differently. Setting x = V^/a 
and y = r]a^, we obtain the equations 



y 



1 

2a; 



In 



1 



1-x 



1] 



(86) 
(87) 



[Wri^{x-1))+Wr{^ix+1))]' 

2 y 

a = -. 

V 

For given a;, we can obtain y from Eq. (j86p [see also Fig. 
[T5] . rj from Eq. ([57)1 and a from Eq. ([55)1 . Varying x in the 

interval ] — 1, 1 [ yields the full curve ric^\a). By symmetry, 
we can restrict ourselves to the interval x € [0,1[. 

The limit a +00 (i.e. y +00) corresponds to 
K a- (i.e. a; ^ 1-). For .t ^ 1", Eqs. (US]) and ((87)) 
reduce to 

y^-l-ln(l-x)^+(x, (89) 



and 



[Wr{^{x-1))+Wr{2^)] 

Using Eqs. (|213p and (|214p . we obtain at leading order 
2 ^ / 1 \ / 1 



(90) 



77 : 



SO that 



1 - -L 

4y 



77i±)(a) 



2 1 



1 

4y 



2 1 



4770^ 



4a' 



+ 



-00 



(91) 



(92) 
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The result ri[^\a) 2 can be understood simply. For 
a — !• +00, the two humps are far away from each other so 

that we expect that the critical temperature rjc^'^ will co- 
incide with the critical temperature of a single Maxwellian 
since they do not "see" each other. Noting that the mass 
of a single hump is M/2, the corresponding critical tem- 
perature is Tc = ^^^^ - ^ leading to 77^=^' = 2. 

On the other hand, the branch r]i^\a) starts at the 
point (a*,?7*), corresponding to ilr = ±K = (i.e. x = 0), 
connecting the branch ric^\a) along which = 0. Now, 



10 



for a; = 0, we have y 
and a* 



1 leading to r/, = = 3.633 



y^Wril) = 0.524. Therefore, the branch ri^c^\a) 
corresponding to f2r = i^* starts at the intersection be- 
tween the branch 77^*^^ (a) corresponding to i?^ = and the 
hyperbole ry = separating the regions where the DF 
has one or two maxima. 
In conclusion: 

• if a > a*, there exists a single critical temperature 
.(±)/ 



riir'{a) determined by Eqs. corresponding to a 

marginal mode {fir = ±K, i^i = 0). Note that the modes 
Qr = +y* and 12^ = — K are degenerate. This degeneracy 
can be raised by a small asymmetry (symmetry breaking) 
in the distribution (see Sec. [6]). 

• if a < a* , there is no marginal mode (Qr — i^* , f^i = 

0). 

The curve r]c^\a) corresponding to the marginal mode 
with pulsation fir = ±14 is plotted in Fig. [181 



5.3 The stability diagram 

The critical temperatures ric{a) corresponding to marginal 
stability determined previously are represented as a func- 
tion of the separation a in Fig. [THl We have also plotted 
the hyperbole 77 = Below this curve, the DF has a 

single maximum at Vb = and above this curve, the DF 
has two maxima at Vq = ±V^, and a minimum at Vq = 0. 
In order to investigate the stability of the solutions in the 
different regions, we can use the Nyquist method. 

For a < a^, there exists two temperatures 77c*'''^^^ 
and 77*;°^'^' at which the DF is marginally stable. For 
r] ~ ?7c°'"^\ the DF has a single maximum at Vb = and 
for 77 = ifc^'^'^^ the DF has a minimum at Vb = and two 
maxima at ±T4- In both cases, the marginal perturbation 
has zero pulsation Qj. = 0- By considering the Nyquist 
curves in this region (see Figs. fT9ll22l) . we find that the 

DF is stable for 77 < r]^^'^^^ and unstable for 77 < 'q^^'^^\ 
there exists one temperature 77c^-' at 



For a > gm 

which the DF is marginally stable. For rj ~ ric^\ the 
DF has two maxima at Vb = ±V^* and one minimum at 
Vb = 0. The marginal perturbation evolves with a pulsa- 
tion Sir = ±y*. By considering the Nyquist curves in this 
region (see Figs. [^^11^5)) . we find that the DF is stable for 

77 < 77c^^ and unstable for 77 > 77^^''. 

For a* < a < qm, there exists three temperatures 
Vc^^''^\ 77c°^^^'' and 77c*'' at which the DF is marginally sta- 



\ .'(y,„J 



Tl=l/a' \ \ 1=n„ (a) 

A \ n =0 
(---) \ \ ' 
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Unstable "-. 
(one mode) \ j-j,_Qj 
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(two modes) 




a=±v. 
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Fig. 18. Stability diagram of the symmetric double-humped 
distribution (|65p . The solid line corresponds to the critical line: 
below this line the DF is stable and above this line the DF is 
unstable. On the left panel (delimited by the solid line and the 
dotted line), there is one mode of instability and on the right 
panel there are two modes of instability. For a* < a < om, 
there exists a "re-entrant" phase as described in the text. The 
dashed line corresponds to the hyperbole rj — below this 

line the DF has one maximum and above this line the DF has 
two maxima and one minimum. The symbols in parenthesis 
give the values of a: or y that parametrize the marginal curves. 
The notations like ( — I — ) give the positions of the er{vext)'s 
in the Nyquist curves as defined in Sec. 12.71 




Fig. 19. Nyquist curve for a < a, and r) < r), 



{0){1) 



< l/a^ 



(specifically a = 0.4 and 77 = 0.5). The DF has only one max- 
imum at Vo = 0. It is stable because the Nyquist curve does 
not encircle the origin. Case (-I-). 



ble. For 77 
maximum at Vq 



„(o)(i) 

'/c 



'iic^^'^\ the DF has a single 



and 77 

0. The marginal perturbation has zero 
pulsation i?^ ~ 0. For 1] — the DF has two maxima 
at Vb = ±T4 and one minimum at Vb = 0. The marginal 
perturbation has a pulsation f2r = ±14 . By considering 
the Nyquist curves in this region (see Figs. I^MHO)) . we 



find that the DF is stable for 7] < rj, 



'/c 



< r] < rj, 



(0)(2) 



stable again for 7/, 



(0)(1) 
(0)(2) 



unstable for 



< 77 < ?7, 



(±) 
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Fig. 20. Nyquist curve for a < a, and t^c^'"'''' < ?7 < 1/a^ < 
77c'^'"'^' (specifically a = 0.4 and ry = 4). The DF has only one 
maximum at Vb = 0. The DF is unstable because the Nyquist 
curve encircles the origin. Case (— ). 



4 - 




-4 - 



-3 -2-10123 
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Fig. 22. Nyquist curve for a < a, and 77 > 77c'^'"'^'' (specifically 
a — 0.4 and i] = 10). The DF has two maxima at Vo = ±V, 
and one minimum at Vo = 0. The DF is unstable because the 
Nyquist curve encircles the origin. Since it rotates twice around 
the origin, this implies that there are N = 2 unstable modes 
{uJr,uJi) with iOi > 0. Case ( — I — ). 




Fig. 21. Nyquist curve for a < a* and < rj < 77'"'"^'' 
(specifically a = 0.4 and r] = 7). The DF has two maxima at 
Vb — ±1/, and one minimum at Vb = 0. The DF is unstable 
because the Nyquist curve encircles the origin. Since it rotates 
only once around the origin, this implies that there is = 1 
unstable mode {uJr,uJi) with uji > 0. Case ( ). 
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Fig. 23. Nyquist curve for a > uai and rj < 1/a^ 
(specifically a = 1 and rj = 0.5). The DF has only one maxi- 
mum at Vb = 0. The DF is stable because the Nyquist curve 
does not encircle the origin. Case (+). 




and unstable again for r] > rjc . This corresponds to a 
re-entrant phase F^. 

In conclusion, considering the stability diagram of Fig. 
[T8l the symmetric double-humped distribution is stable 
below the solid line and unstable above it. When the solu- 
tion is unstable, we expect that the system will evolve 
through violent relaxation towards a QSS that corre- 
sponds, in case of efficient mixing, to a Lynden-Bell sta- 
tistical equilibrium state. 



Interestingly, a re-entrant phase has also been observed in 
the stability diagram of the Lynden-Bell (or Fermi-Dirac) dis- 
tributions for the HMF model |55l58j . 




0.5 1 1.5 

e 

r 

Fig. 24. Nyquist curve for a > um and < rj < jji*' 

(specifically a = 1 and rj — 2). The DF has two maxima at 
Vb = ±K, and one minimum at Vb = 0. The DF is stable 
because the Nyquist curve does not encircle the origin. Case 
(+ + +)■ 
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Fig. 25. Nyquist curve for a > om and r] > 77c*' > 
(specifically a = 1 and 77 = 5). The DP has two maxima at 
Vb — i:V, and one minimum at Vb = 0. The DF is unstable 
because the Nyquist curve encircles the origin. Since it rotates 
twice around the origin, this implies that there are N = 2 
unstable modes {u)r,u!i) with cui > 0. Case ( — \ — ). 




Fig. 28. Nyquist curve for a, < a < aM and t]!"'"^'' < tj < 
1/a^ < 77^=^' (specifically a = 0.535 and 77 = 3.4). The DF 
has only one maximum at Vb = 0. The DF is stable because 
the Nyquist curve does not encircle the origin. Case (+). Note 
that when 77 — > , the tangent at {€r{uJr = 0), 0) is no more 
vertical since /"(O) = (see remark at the end of Sec. 12. 5p . 
This corresponds precisely the transition between one maxi- 
mum and two maxima and a minimum. 



uj- 




Fig. 26. Nyquist curve for a, < a < am and 77 < 77'"'"^'' < l/a^ 
(specifically a — 0.535 and 77 — 0.5). The DF has only one 
maximum at Vb = 0. The DF is stable because the Nyquist 
curve does not encircle the origin. Case (+). 



u" 




Fig. 27. Nyquist curve for a, < a < om and 7;'"''^^ < 77 < 
^(0)(2) ^ ^^^2 (specifically a = 0.535 and 77 = 2.5). The DF 
has only one maximum at Vb = 0. The DF is unstable because 
the Nyquist curve encircles the origin once. There is iV = 1 
mode of instability. Case (— ). 
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Fig. 29. Nyquist curve for a* < a < um and 77c'^'''^' < l/a^ < 
rj < rji"^^ (specifically a = 0.535 and 77 = 3.55). The DF has 
two maxima at Vb = ±V, and one minimum at Vb = 0. The 
DF is stable because the Nyquist curve does not encircle the 
origin. Case (+ + +). 



6 The asymmetric double-humped 
distribution 

6.1 Determination of the extrema 

Wc consider an asymmetric double-humped distribution 
of the form 



/(«) 



+ Ae~ 



(93) 



where T = is the temperature of the MaxweUians and 
A is the asymmetry parameter (we assume here that A > 
1). This distribution is plotted in Fig. [511 The symmetric 
case is recovered for A = 1. The average velocity is {v) ~ 
— [{A — 1)/{A + l)]va and the kinetic temperature Tkin = 
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Fig. 30. Nyquist curve for a, < a < aM and i] > 7?i*' > 1/a^ 
(specifically a = 0.535 and 77 = 5). The DP has two maxima 
at Vo = ±14 and one minimum at Vb = 0. The DF is unstable 
because the Nyquist curve encircles the origin twice. There are 
N = 2 modes of instability. Case ( — I — ). 
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Fig. 31. Asymmetric double-humped distribution made of two 
Mawellians with separation a and asymmetry Z\ > 1. If > 
y*/a?, the DF has one global maximum at Vb = VI < 0, a 
minimum at Vb = V^ > and a local maximum at Vo = V+ > 
0. If < y*/a? , the DF has only one maximum at Vb = V'- < 0. 



{{v - {v)f) = T+ [4Z\/(Z\ + l)^]vl. The velocities vq at 
which the distribution function f{v) is extremum satisfy 
f'{vo) = 0. They are determined by the equation 



(94) 



Introducing the dimensionlcss temperature (|37[) . the di- 
mensionless velocity psp and the dimensionless separation 
(p7)) . Eq. dnil) can be rewritten 



2aVn 



■In 



Vo 



a-Vn 



MA) 

2aVn ' 



(95) 



For a given value of inverse temperature 77, separation a 
and asymmetry A, this equation determines the velocities 
Vo where /(V) is extremum. We note that Vo g] — a, +a[- 
It is convenient to introduce the variables 



Vo/a, 



y = m 



(96) 



For fixed a, the parameter y plays the role of the inverse 
temperature and the parameter x plays the role of the 
velocity. Then, we have to study the function 



1 

22; 



■In 



1 



1 -a; 



ln(Z\) 

2x ' 



for X s] — 1, +1[. This function is plotted in Fig. 
the following properties 



y{x) - 



4 Ml 



1 



■ln(l + j; 

\n{A) 
2x ' 



(x^l-), 
{x ^ 0). 



(97) 
It has 

(98) 
(99) 
(100) 



Considering the negative velocities a; < 0, we note that 
y > iS X < xq with 



Xo 



A~l 

aTi' 



(101) 
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\ (x^.y^) if A < A, 
X^^^^^^^x^.y^) if A>^A^^ 
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Fig. 32. The function y{x) for the asymmetric double-humped 
distribution. 



Considering the positive velocities a; > 0, we note that the 
curve y{x) is minimum at x = x^ where x^, is solution of 



2x^ 



1 



In 



1 - a;. 



In(zi). 



(102) 



This fimction is represented in Fig. [331 We note that Eq. 
(|102p has a unique solution a;* for each value of Z\ > 1. 
Therefore, the function y{x) has a single minimum at a; = 
x^ . The value of this minimum is 



1 



1 



> 1. 



(103) 



The extrema of the distribution f{V) can be determined 
from the study of the function ((97l) . If j/ < y*, i.e. 
V < y*l'^ 1 the distribution /(V^) has a single maximum 
at Vb = T^- < 0. If 2/ > y*, i.e. r/ > y^jc? , the distribution 
/(F) has two maxima at Vb = F_ < and Vb = V+ > 
and one minimum at Vq = V^ > 0. These different values 
arc given by Vb = ay~^(r\c?\ 

In conclusion, for a given asymmetry A and separation 

a: 
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Fig. 33. Evolution of x,, y,, Xc and yc as a function of A. The 
curves intersect each other at A, — 3.3105. 



• if 77 > ?/*/a^, the distribution function f{V) has a 
global maximum at Vq = < 0, a local maximum at 
Vb = V+ > and one minimum at Vq = Vp > 0. 

• if 77 < ytfja}^ the distribution function /(V) has only 
one maximum at Vq = 1^- < 0. 

Remark: We can note the formal analogy with the 
mean field theory (ferromagnetic transition) of the one 
dimensional Ising model with magnetic field where the 
asymmetry A plays the role of the magnetic field (com- 
pare Eq. jMl) with Eq. (14) of [86]). 



6.2 The condition of marginal stability 

The dielectric function associated to the asymmetric 
double-humped distribution is 

<n) = 1 - [WiVvi^ - a)) + AW{^i{n + a))] , 

(104) 

where W{z) is defined in Eq. (|40)) . When i7j = 0, the real 
and imaginary parts of the dielectric function e(i7r) = 
er{fir) + i^i{^r) Can be written 

er[Qr) = 1 - -Y^iWriVvi^r - a)) 

+AWriy/v{fir + a))], (105) 



-[W,{^{Qr - a)) 



l + A' 

+Z\W,(V^(a + a))], (106) 

where Wr{z) and Wi{z) are defined in Eqs. ((44|) - ([45|) where 
z is here a real number. The condition of marginal stabil- 
ity corresponds to e,.(i7,.) = ei(J7r) = 0. The condition 
ei{Qr) = is equivalent to 

n^r) = 0. (107) 

The condition er{fir) — leads to 

1 - [Wr{^^{Qr - a)) + AWr{^{[2r + «))] = 0. 



(108) 



According to Eq. (jl07p . the real pulsation ilr is equal to 
a velocity Vg at which the distribution (j93|) is extremum. 
The second equation (|108p determines the value(s) ric{a) 
of the temperature at which the distribution is marginally 
stable. Therefore, we have to solve 



l + A 



[Wr{^/T]{V,, - a)) + AWr{^i{V,, + a))] = 0, 



where Vq is given by 
1 



'7 



2aVb 



In 



In(^) 
2aVo ■ 



(109) 



(110) 



Eliminating Vq between these two expressions yields the 
critical temperature(s) ?/c(a) as a function of a. However, 
it is easier to proceed differently. Setting a; = Vo/a and 
y = ija? , we obtain the equations 



1 l + x 
l + A 



In(zi) 
2x '■ 



(111) 



(112) 



(113) 



[Wr{^{x - 1)) + AWr{^{x + 1))] ' 

2 y 

a = —. 
V 

For given x, we can obtain y from Eq. (|llip [see also Fig. 
[32] . 77 from Eq. pT2]l and a from Eq. (fTTS]) . Varying a; in 
the interval ] — 1, 1[ yields the full curve ijda). We have 
three types of solutions. For x g] — l,Xo], we obtain a 
branch 77c ^ (a) where the pulsation of the marginal mode 
is negative: fir = V- < corresponding to the global 
maximum of f{v). For x g]0, x*], we obtain a branch 

7?i^'' (a) where the pulsation of the marginal mode is posi- 
tive: f2r = Vp > corresponding to the minimum of f{v). 

For x e [x,, 1[, we obtain a branch 77c^^(a) where the pul- 
sation of the marginal mode is positive: f2r = V+ > 
corresponding to the local maximum of f{v). This leads 
to the curves reported in Figs. [34l [35| and [36| for three 
values of the asymmetry factor A. 

The limit a — > +00 (i.e. y +00) corresponds to 
Vq (i.e. X 1^) or Vb ^ — a+ (i.e. x —1^)- For 

X ^ 1", Eqs. (fTTT]l and (fTT^ reduce to 



and 



y iln(l-a;) +00, 



l + A 



(114) 
(115) 



[Wriy^ix ~ 1)) + AWri2^)]' 

Using Eqs. (j213p and (|214p . we obtain at leading order 



l±|.a..,(..|).a..,(. 



A 



so that 



77i+)(a) ^l + A + 



A 

4^ 



(116) 



(117) 
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Fig. 34. Stability diagram of the asymmetric double-humped 
distribution for A < Ac < Zi, (specifically A = 1.02) showing 
a re-entrant phase in continuity with the symmetric case. Since 
A < Z\*, there exists three marginal branches rji \a), rji^\a) 
and 77c''''' (a). 
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Fig. 35. Stability diagram of the asymmetric double-humped 
distribution for Ac < A < A* (specifically A = 1.5) showing 
the disappearance of the re-entrant phase. 
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Fig. 36. Stability diagram of the asymmetric double-humped 
distribution for A > A* (specifically A = 10). Since A > A», 
there exists only two marginal branches r)\r\a) and ri^\a). 



For a: ^ Eqs. (fTTT|) and (fTT^ reduce to 

1 



y ' 



■ ln(l -I- x) ^ +00, 



and 



(118) 
(119) 



l + A 

Using Eqs. (|213p and (|214p . we obtain at leading order 



l + A l + Af_^, 1 



A-i- A 



AAy 



l + A 
A 



I 



AA'qa'^J ' 
(120) 



so that 



l + A 
A 



AAa? 



(a^+oo). (121) 



The limit a — > with rj finite corresponds to y — > 0, 
i.e. a; ^ xo. We find from Eq. (fm]) that -q l/Wr{0) = I 
so that 

#)(a)^l (a^O). (122) 



This retm'ns the critical temperature (|46l) associated with 
the Maxwellian distribution ([34]) corresponding to a = 0. 
The limit a — > with ?/ — !■ +oo corresponds to 



(a^O), 



(123) 



with Uc finite {P ^ p ii A < A^, and P = + if Z\ > zi* as 
will become clear below). Since rj +oo, this constant yc 
is determined by 



1 , fl + Xc 

Vc = - — hi 



I — Xc 



2xc ' 



(124) 



Wr{^c{Xc - 1)) + AWr{^c{x, + I)) = 0. (125) 

Note that there is no physical solution to Eqs. ()llip - (jll3p 
when < a; < a;c (?7 would be negative) so that the branch 
r]i^\a) starts at (0,+oo) corresponding to a; = Xc- The 
evolution of Xc and j/c with A is studied in Fig. [33] and 
is compared with the evolution of a:* and j/*. The curves 
intersect each other at A^, = 3.3105. For A < A^,, Xc < 
X* so that the stability diagram displays three marginal 

branches r/i \a), r{f\a) and •q"c^\a) (see Fig. [34]). The 
branches 77^^^ (a) and jyc^-* (a) connect each other at (a* , 77* ) 
corresponding to a; = a;*. At that point they touch the 
hyperbole 77 = y*/a^ separating distributions with one or 

two maxima. Since j/c > J/*, the curve 'r]c^\a) for a ^ 
is above the hyperbole 77 = y^./a'^. For A = A^,, Xc = 
a;* so that the branch 77^^'' (a) is rejected at infinity and 
only the branches 77c \a) and 77c^^(a) remain. Since t/c = 

77*, the branch 770^'* (o) coincides with the hyperbole 77 = 
77*/a^ for a ^ 0. For A > A^, Xc > x^ so that the phase 
diagram displays only two marginal branches 77c ^ (a) and 

77c^-'(a) (see Fig. [351) • The curve 77c'''^(a) is stictly above 
the hyperbole rj ~ y^ja? without any intersection. 
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6.3 The stability diagram 

The critical temperatures r]c{a) corresponding to marginal 
stability determined previously are represented as a func- 
tion of the separation a in Figs. [34l[35l and[36lfor three val- 
ues of the asymmetry factor A. We recall that Tyi"' (a) cor- 
responds to the temperature associated with a marginal 
mode with pulsation f2r = V- < 0, ■q"c\a) corresponds 
to the temperature associated with a marginal mode with 
pulsation fir = Vp > Q and 7{c^\a) corresponds to the 
temperature associated with a marginal mode with pul- 
sation fir = V+ > 0. We have also plotted the hyperbole 
77 = y^,/a^ ■ Below this curve, the DF has a single maximum 
at y_ < and above this curve, the DF has two maxima 
at y_ < and V+ > and a minimum at > 0. In order 
to investigate the stability of the solutions in the different 
regions of the parameter space, we have used the Nyquist 
criterion. There exists a critical value Ac = 1.09 of the 
asymmetry parameter. For A < Zic, the stability diagram 
displays a re-entrant phase (in continuity with the case 
Zi = 1) but for A > Ac the re-entrant phase disappears. 
The description of the stability diagram is similar to the 
one given in Sec. 15.31 and the different possible cases can 
be understood directly from the reading of Figs. [34l [35] 
and 1361 The best way is to fix a and progressively increase 
the value of t]. Below the dashed line, the distribution has 
only one maximum at Vq — V- so the Nyquist curve has 
one intersection with the x-axis (in addition to the limit 
point (1,0)). For i] 0, we find that er(VL) > so the 
system is stable. As we increase i] and pass above the 
dashed line, the distribution has two maxima at Vq = V- 
and Vq = V+ and one minimum at Vb = At each in- 
tersection with a marginal line, one of the values er{V-), 
er{Vp) or e, (V+) changes sign. We have indicated by sym- 
bols like ( h +) the repsective signs of er{V-), er{Vp) 

and er{V+). We can then easily draw by hands the corre- 
sponding Nyquist curve. Therefore, it is not necessary to 
show all the possibilities and we have only indicated a few 
representative cases in Figs. I37ll39l for illustration. 



7 The repulsive HMF model 
7.1 General results 

We shall now consider the HMF model with a repulsive 
interaction between particles. It is described by the Hamil- 
tonian 

1=1 i<j 

with fc > 0. In this paper, we shall investigate the linear 
dynamical stabilityEj of stationary solutions of the Vlasov 



The nonlinear regime of the repulsive HMF model leading 
to a bicluster (for low energies) has been studied in [87] . 



4 



2 



-2 



-4 





1 


' 1 ' 


y v^ 








y p 






1 


, 1 , 



Fig. 37. For A < Ac < At: Nyquist curve for a < a* and 
j/,/a^ < 77 < rii^^ (specifically A = 1.02, a = 0.4 and rj = 8). 
The DF has two maxima at V- and V+ and one minimum at 
Vp. The DF is unstable because the Nyquist curve encircles 
the origin once. There is TV = 1 unstable mode. Case ( ). 
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Fig. 38. For Ac < A < A*: Nyquist curve for yt/a^ < rj < 
(specifically A = 1.5, a = 1 and 77 = 1.8). The DF has 
two maxima at V- and V+ and one minimum at Vp. The DF is 
stable because the Nyquist curve does not encircle the origin. 
Case (+ + -f). 
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Fig. 39. For Ac < A < A,: Nyquist curve for y,/o? < 77^ ' < 
ri < rji'^^ (specifically A = 1.5, a = 1 and 77 = 2.5). The DF 
has two maxima at V- and V+ and one minimum at Vp. The 
DF is unstable because the Nyquist curve encircles the origin 
once. There is A'^ = 1 unstable mode. Case ( — h +). 



P.H. Chavanis and L. Delfini: Application of the Nyquist method to the HMF model 



21 




Fig. 40. For Ac < A < A^: Nyquist curve for 7] > ric ' ' > 
y»/a^ (specifically A — 1.5, a — 1 and 77 = 6). The DF has two 
maxima at V- and V+ and one minimum at Vp. The DF is 
unstable because the Nyquist curve encircles the origin twice. 
There are N = 2 unstable modes. Case ( — \ — ). 



equation of the form / = f{v). In the repulsive case, the 
dispersion relation becomes 



k f f'iv) 
e{uj) = 1-- I ^!-^dv ^ 0. 



(127) 



When uji ~ 0, the real and imaginary parts of the 
dielectric function €{10^.) = Sri^r) + i^i{^r) ar^ 



k f 



+00 (•/ 



V — UJr 



(128) 



(129) 



To apply the Nyquist method, we have to plot the curve 
(eriiOr), eiiuJr)) parameterized by ujr going from — cx) to 
+00. Let us consider the asymptotic behavior for lu^. — > 
±00. Since f{v) tends to zero for v ±00, we conclude 
that ei{ujr) — > for Ur —>■ ±00 and that ei{u>r) < for 
LOr — * —00 while £i{iOr) > for LUr — *■ +00. On the other 
hand, for ujr — > ±00, we obtain at leading order 

k p 
2^ 



er(Wr) 



1 P / 



±00). 



(130) 



From these results, we conclude that the behavior of the 
curve close to the point (1, 0) is the one represented in Fig. 
1411 Let Vext be the velocity corresponding to an extremum 
of the distribution (f'{vext) = 0). Then, we have 



2 ./ „ iv - 



/(«) 



dv > 1. (131) 



If VMax denotes the velocity corresponding to the global 
maximum of the distribution, we clearly have 



(132) 



Finally, for a repulsive interaction, we note that if the zero 
of f'{v) corresponds to a maximum (resp. minimum) of / 
then the hodograph crosses the real axis upward (resp. 
downward) . 
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Fig. 41. Nyquist curve for the Maxwell and the Tsallis distri- 
bution with 71 < —1 (specifically n = —2). 



uj" 




Fig. 42. Nyquist curve for the Tsallis distribution with n > 
5/2 (specifically n — 5). 



7.1.1 Single-humped distributions 

Let us assume that the distribution f{v) has a single max- 
imum at u = vq (so that /'(wo) = 0) and tends to zero for 
V ±00. Then, the Nyquist curve cuts the x-axis (6^(0;^) 
vanishes) at the limit point (1,0) when ujr ±oo and 
at the point (er(i'o),0) when ujr = vq. Due to its behavior 
close to the limit point (1,0), the fact that it rotates in the 
counterclockwise sense, and the property that er(wo) > 1 
according to Eq. (|132[) , the Nyquist curve must necessarily 
behave like in Fig.[4TJ Therefore, the Nyquist curve starts 
on the real axis at er{LOr) — 1 for ujr —00, then going in 
counterclockwise sense it crosses the real axis at the point 
^rivo) > 1 and returns on the real axis at er(ct'r) = 1 for 
LOr +00. Therefore, it cannot encircle the origin. Ac- 
cording to the Nyquist criterion exposed in Sec. 12.41 we 
conclude that a single-humped distribution is always lin- 
early stable. For illustration, the Nyquist curves for the 
Maxwell and the Tsallis distributions are represented in 
Figs.SIEl 

For the water-bag model (ri = 1/2), we can obtain an 
analytical expression of the dielectric function in the form 



(133) 
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Fig. 43. Nyquist curve for the Tsallis distribution with 3/2 < 
n < 5/2 (specifically n — 2). 




Fig. 44. Nyquist curve for the Tsallis distribution with 1/2 < 
n < 3/2 showing the difi'erent cases (specifically n = 3/2, n — 
1.2, n = 1, n = 0.7 and n = 1/2). 



where W'^^/^^z) is given by Eq. (|6ip . The condition 
e{uj) = determines the complex pulsation. We have 



kM 
An 



3T. 



(134) 



The pulsation is real so the mode is purely oscillating with- 
out being damped. For T = 0, we find that uj = ujp where 
ujp is the proper pulsation defined in Eq. (|25ip . Note, in 
passing, that the dimensionless pulsation can be written 
i? = Lu/ujp. On the other hand, for f2i = 0, we get 



Nyquist curves starts at (1,0), progresses in the counter- 
clockwise sense and crosses the a;- axis at tri^uax) > 1, 
then at er{vmin) and er{vmax)- We can convince ourselves 
by making drawings of the following results. If 

(+ + +): erivMax) > 0, eriVmin) > 0, eriVmax) > 0, 

(H ): er(vMax) > 0, CriVmin) < 0, CriVn^ax) < 0, 

(+ + -): erivMax) > 0, er(Wmm) > 0, er{Vmax) < 0, 

the Nyquist curve does not encircle the origin so the sys- 
tem is stable. If 

(H h): er{vMax) > 0, CriVmin) < 0, 6^ ( )>o, 

the Nyquist curve rotates one time around the origin so 
that there is one mode of instability. Since er{vMax) > 
there is no mode of marginal stability with tUr = VMax- 

Cases (+ + +), (H h) and (H ) are observed in Sec. 

17.31 for an asymmetric double-humped distribution made 
of two Maxwellians. 

If the double-humped distribution is symmetric with 
respect to the origin with two maxima at ±t;* and a min- 
imum at w = 0, we get the same results as above with the 
additional properties er{vMax) = ^rivmax) = er(i'*) > 1 
and er{v„iin) = er(0). We have only two cases (+ + +) 
and (H h). They are observed in Sec. l7.2l for a symmet- 
ric double-humped distribution made of two Maxwellians. 
Since er{v*) > 0, there is no mode of marginal stability 
with cj,. = zbw* . 



7.2 The symmetric double-humped distribution 

Let us now consider the symmetric double-humped distri- 
bution (j65p . In the repulsive case, the dielectric function 
is 

e{f2) = 1 + 1 [W{^{Q - a)) + W{^{Q + a))] . (136) 

The condition of marginal stability corresponds to e(I?) = 
and LOi = 0. The condition ei{LUr) = is equivalent 
to f'{LOr) = so that the real pulsations correspond to 
the velocities where the distribution is extremum: lo,. ~ 
vq. Then, the temperatures at which the distribution is 
marginally stable are obtained by solving er{LOr — vq) = 0. 
Proceeding as in Sec. [5] and introducing the parameters 
X = Vo/a and y = ija^, we find that the equations deter- 
mining the critical temperatures rjda) are given by 



^lnfi±^ 
2x \l~x 



(137) 



1 + 



3' I- 



e^{nr) = 0. (135) 



Therefore, the Nyquist curve is made of two segments 
oo, 1] and [1 -I- 77/3, -|-oo[ as represented in Fig. O 



7.1.2 Double-humped distributions 

Let us consider a double-humped distribution with a 
global maximum at Wj\/aa: , a- minimum at Vmin and a local 
maximum at Vraax- We assume VMax < Vmin < Vmax- The 



V 



(138) 



(139) 



[Wr{^ix~l))+Wr{^ix+1))]' 

2 y 

a = —. 

V 

We note that only the sign in Eq. (|138p changes with re- 
spect to the study of the attractive case, so we can readily 
adapt the results of Sec.[5]to the present situation by sim- 
ply reverting the sign. For i?,. = ±V^,, corresponding to 
X ^ 0, there is no physical solution to Eqs. (|137p - (|139p 
with positive temperature 77 > 0. Therefore, in the repul- 
sive case, there is no marginal mode with non zero pul- 
sation for the symmetric double-humped distribution (in 
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Fig. 45. Stability diagram of the symmetric double-humped 
distribution for the repulsive HMF model. There is a re-entrant 
phase for 1 < a < am- 
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Fig. 46. Nyquist curve for a < 1 and r) < Ijc? < 77'"'"^'' (a) 
(specifically a = 0.5, ry = 2). The system is stable. Case {+). 



agreement with the general discussion of Sec. I7.1.2[) . We 
now consider the marginal mode with i?^ = 0. This corre- 
sponds to the "degenerate" solution a; = (for any y) for 
which Eqs. p37)) - p39)) reduce to 



= -yWriVy)- 



(140) 



(141) 



According to Fig. I16[ physical solutions exist only for y > 
Umax = -^c- We note that the range of parameters that 
was forbidden in the attractive case is now allowed in the 
repulsive case and vice versa. We thus consider the range 
y G [ymaxi +00 [. The function p41|) has a maximum a„i = 
1.282 at y = ?/,„ = 7.642 and for y +00 it tends to 
a = 1. Therefore, for a > am, Eqs. (|140p - (|14ip have no 
solution. For a = a,n, Eq. (|14ip has one solution y = 
determining, through Eq. (|140|) . a temperature r]i°^^"^^ 
4.647. For 1 < a < a™, Eq. (|141[) has two solutions y = 

yifi- Then, Eq. ()140p determines two temperatures 77c"'"^'' 

and ?yc°''^^'* as illustrated in Fig. [T71 Finally, for < a < 1, 
Eq. (|14ip has one solution y ^ yi determining, through 

Eq. ()140|1 . a temperature r?i°^^^^ 

The limit y — > ymax corresponds to a ^ 0. Since ~ 
-ymaxWriy/y) and r; -l/Wr{,/y), we get 



V, 



(0)(1) 



(«) - 7l> 



(a-.0). 



(142) 



This result is to be expected since, for a = 0, the distri- 
bution 165]) reduces to the Maxwellian that is stable for a 
repulsive interaction (hence 77c = +00). 

Using Eq. (j213p . the limit y s- +00 corresponds to 
a ^ 1+. Then, we get 



V, 



(0)(2) 



(a) 



2(a- 1)' 



(a^l+). 



(143) 



The critical temperature(s) 77c (a) corresponding to 
marginal stability determined previously are represented 
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Fig. 47. Nyquist curve for a < 1 and < 77 < 77c"'"^'(a) 

(specifically a — 0.5, 77 = 5). The system is stable. Case (+H — h). 



as a function of the separation a in Fig. 351 We have also 
plotted the hyperbole 77 — 1/a^. Below this curve, the dis- 
tribution has a single maximum at Vq = and above this 
curve, the distribution has two maxima at Vb = ±K and 
a minimum at Vq = 0. In order to investigate the stability 
of the solutions in the different regions, we have used the 
Nyquist criterion. 

For a < 1, there exists one temperature 77c'' ''^^^ at which 
the distribution is marginally stable. For 77 = 77i'^^'"'^'', the 
distribution has a minimum at Vq = and two maxima at 
±T4- The marginal perturbation has zero pulsation fi^ = 
0. By considering the Nyquist curves in this region (see 

Figs. HSlHSl) . we find that the DF is stable for 77 < 77^"'*^^ 
and unstable for 77 > rji'^'^^^^ . 

For a > ttrrn there is no marginal mode and the dis- 
tribution is always stable (see Figs. H51l5(I)) . This result 
is to be expected because, for a — > -l-oo, the two humps 
do not "see" each other and behave as two independent 
single-humps distributions that are stable in the repulsive 
case. 



(0)(1) 



For 1 < a < a,„, there exists two temperatures 77, 
and 77c'^''^'' at which the distribution is marginally sta- 
ble. These two branches merge at the point («„, 77^*^^ '■'"■'). 
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Fig. 48. Nyquist curve for a < 1 and 77 > Vc'^ '{a) > Fig. 51. Nyquist curve for 1 < a < am and rj < 1/a (specifi- 

(specifically a = 0.5, rj = 15). The system is unstable. Case cally a = 1.2, rj = 0.5). The system is stable. Case (+). 
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Fig. 52. Nyquist curve for 1 < a < Um and < rj < 

Fig. 49. Nyquist curve for a > and rj < Ijc? (specifically rjf"'^^\a) (specifically a = 1.2, rj = 3). The system is stable. 
a = lA, rj = 0.4). The system is stable. Case (+). Case (+ + +). 
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Fig. 50. Nyquist curve for a > am and rj > 1/a^ (specifically 
a = 1.4, rj = 1). The system is stable. Case (+ + +). 



For 77 = jyc"''^^ and j] — rjc^^^^\ the DF has a minimum 
at Vo = and two maxima at ±T4. The marginal per- 
turbation has zero pulsation f2r — 0. By considering the 
Nyquist curves in this region (see Figs. [CTIIS^ . we find 

that the distribution is stable for 77 < Tji^'''''^\ unstable for 
j^(o)(i) ^ ^ 77c'^'""^-' and stable again for 77 > rjc^^'"^\ This 
corresponds to a re-entrant phase. 

7.3 The asymmetric double-humped distribution 

The dielectric function associated to the asymmetric 
double-humped distribution ([M)) in the repulsive case is 

<n) = 1 + [WiVvi^ - a)) + AWiVTj{Q + a))] . 

(144) 

Proceeding as in Sec. [6] and introducing the parameters 
X = Vq/u and y = rjo? , we find that the equations deter- 
mining the critical temperatures rjc{a) are given by 

1 /l + a;\ ln(Z\) 
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Fig. 53. Nyquist curve for 1 < a < a,n and 1/a < ?7c (a) < 
»7 < 77c°'"'^'' (a) (specifically a = 1.2, r] — 7). The system is 
unstable. Case (H h). 




Fig. 54. Nyquist curve for 1 < a < am and t] > rjc (a) > 
(specifically a — 1.2, rj — 12). The system is stable. Case 
(+ + +)■ 
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Fig. 55. Stability diagram for the asymmetric double-humped 
distribution with A < At < Ai^^ (specifically A = 2). We have 
also plotted the hyperbole 77 — ytja^ ■ Below this curve, the 
distribution has a single maximum at > and above this 
curve, the distribution has a global maximum at Vb = VI < 0, 
a local maximum at Vb = V+ > and a miinmum at > 0. 
Since A < Zi* , there exists only one marginal branch r]^^ which 
is always stricly above the hyperbole r\ — y»/a? . 
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Fig. 56. Stability diagram for the asymmetric double-humped 
distribution with A, < A < A^c^ (specifically A = 20). Since 



A > At,, there exists two marginal branch ini^J'^ and r))^' that 
connect each other at the point of contact with the hyperbole 
7? = y,/a^. 



„(+) 



= ^. (147) 
V 

Equations (|145|) and ()147|) determine the extrema of the 
distribution and Eq. p46p determines the temperatures 
corresponding to the modes of marginal stabihty. As in 
Sec. [6l the curve r]c{a) can be obtained by varying x 
between —1 and -1-1. In the repulsive case, there ex- 
ists physical solutions with positive temperature only for 
< a; < a;c. For A < A^ , we get only one marginal branch 
rjc^^ (a) corresponding to the mode Hr = Vp > (see Fig. 
[55|) . For A > we get two marginal branches 'qi'\a) 
and 77c^^ (a) corresponding to the modes i7,. = > and 
i7r = ^+ > (sec Fig. [56|) . They connect each other at 
(a*, 77*) corresponding to a; = x*. At that point they touch 
the hyperbole 77 = y^l c? separating distributions with one 
or two maxima. 



For a; ^ 0+, using Eq. ((2T3)) . we find that a \. 
Then, proceeding carefully, we get 

For a; — *■ a;c, we find that a ^ and 

77f)(a)^^, (a^O). (149) 

The stability diagrams corresponding to the asymmet- 
ric double-humped distribution with A =- 2 and = 20 
are represented in Figs. [55l and 1561 In continuity with the 
symmetric case Z\ = 1, they display a re-entrant phase. 
Some representative Nyquist curves are represented in 
Figs. [57115^ In Fig. [SDl we determine how the stability 
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Fig. 58. For Z\ < Z\* < A'^P : Nyquist curve for j/,/a^ < 
ri'f\a) < ri (specifically A = 2, a = 0.9 and ?? = 10). The 
system is unstable. Case (H h). 
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Fig. 60. Stability diagram for different values of the asymme- 
try A showing the appearance of a double re-entrant phase at 
Zic^' = 25.6 and their disappearance above Zic^' = 43.9. 



diagram changes with the asymmetry A. We find that a 
double re-entrant phase appears for A'^^ — 25.6 and that 

(2) 

they both disappear above Ac ~ 43.9. In fact, the critical 
value Ai^^ corresponds to tlie case where the asymptote 
passes from a = l"'"toa=l~. Therefore, it corresponds 
to the vanishing of the term in parenthesis in Eq. p48p . 
A numerical solution of 3 + In Z\(l - Z\)/(l + Zi) =0 then 

gives a'-c^ = 25.6268.... 



8.1 Linear stability for the barotropic Euler equation 

We consider the dynamics of a gas governed by the Euler 
equations 



dp ^ 
— + V 
dt 



(pu) - 0, (150) 



du 



V)u 



^-Vp-pV$, (151) 



8 General results for an arbitrary potential of 
interaction 

In this section, we show that the previous results obtained 
for the attractive or repulsive HMF models can have appli- 
cations for more general potentials of interaction. Wc also 
make the connexion between linear and formal nonlinear 
stability, and between dynamical stability with respect to 
the the Euler and Vlasov equations. 



<?(r,t) = J ui\r~r'\)pir',t)dr', (152) 

where u(|r— r'|) is an arbitrary binary potential of interac- 
tion. We assume that the gas is described by a barotropic 
equation of state p ~ pip)- Clearly, a spatially homoge- 
neous distribution p = <P = Cst. with u = is a station- 
ary solution of Eq. p50|l - (|152p provided that <P ^ pU with 
U ~ J u(x)rfx (for the gravitational potential, this last 
condition is not realized and wc have to make the Jeans 
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swindle [21]). Wc want to investigate the linear dynam- 
ical stability of this spatially homogeneous distribution. 
Linearizing the Euler equations p50[) - (|152p around this 
steady state and decomposing the perturbations in nor- 



mal modes Sp^^j ~ iJ^kw 
the dispersion relation [29] 

iu^ = elk' 



Sui^ui ''^ e^^^''^ we obtain 



(153) 



where we have introduced the velocity of sound = p'{p)- 
We note that the pulsation lu is either real (w^ > 0) or 
purely imaginary (w^ < 0). As a result, the system is lin- 
early stable with respect to a perturbation with wavenum- 
ber k if 

cl + {2TTfu{k)p > 0, (154) 

and linearly unstable otherwise (when uj"^ < 0, the mode 
uj = +zv^Lj^ grows exponentially rapidly). 

For repulsive potentials satisfying u{k) > 0, the ho- 
mogeneous distribution is always stable. For Coulombian 
plasmas in d = 3, using Eq. p40p . the dispersion relation 
can be written [721 : 



2; 2 

c„k , 



where ojp is the plasma pulsation Eq. (|24ip . 



(155) 
For the re- 



pulsive HMF model, using m„ = ■^Sn,±i, 
relation is [29] : 

2 2 2, ^-^'^ 2 r 

UJ ^ c^n + - — n dn.±i. 
47r 



the dispersion 



(156) 



±Cs?7.. 

^ c?„ where 



The modes n 7^ ±1 oscillate with a pulsation 
The pulsation of the modes n — ±.\\s, uP' = lo . 
ujp is the proper pulsation (|25ip . 

For attractive potentials satisfying u{k) < 0, the ho- 
mogeneous distribution is linearly stable if 

E {27Tfp\u{k)U 



(157) 



and linearly unstable if < {cDcrit- In this last case, 
the unstable wavelengths are determined by the converse 
of inequality (|154p . For the gravitational interaction in 
d = 3, using Eq. (|246p . the dispersion relation can be 
written [21] : 

(158) 
00) for suffi- 

(159) 



^2 = c^fc2 - AttGp- 

The system is always unstable ((c^), 
ciently small wavenumbers 



1/2 



where kj is the Jeans wavenumber for a barotropic gas. 
The growth rate is maximum for fc 0. For the attractive 
HMF model, using m„ = —^Sn ±1, the dispersion relation 
is (Ml: 



2 2 



kM 
An 



■n 6n. 



±1- 



(160) 



The modes n ^ ±1 oscillate with a pulsation lo = ±Csn. 
The complex pulsation of the modes n = ±1 is = 
cl - The system is stable for > and unstable 
for cl < M£. 



8.2 Formal stability for the barotropic Euler equation 

The barotropic Euler equation conserve the mass M = 
J pdr and the energy functional [21] : 



W[p, u] = J pJ'P^ dp'dv + lJp<Pdr 



p— dr. 



(161) 



Therefore, a minimum of the energy functional at fixed 
mass determines a steady state of the Euler equation that 
is nonlinearly dynamically stable [3- The qualitative idea 
is the following. If the system is an energy minimum, and 
since the energy is conserved by the Euler equation, then 
it cannot evolve away from that point. Therefore, it is 
dynamically stable. Alternatively, if the flow is a saddle 
point of energy, then it can evolve away from that point 
by following an iso-energy line. In that case, it is not dy- 
namically stable. Here, we have stability in the sense of 
Lyapunov. This means that the size of the perturbation 
is bounded by the size of the initial perturbation for all 
times. We are led therefore to considering the minimiza- 
tion problem 



min {VV[p, u] 



M[p] = A/}. 



(162) 



In this paper we shall consider small perturbations and 
determine conditions under which the system is a local 
minimum of energy {SW = and S'^W positive definite). 
Thus, we restrict ourselves to formal nonlinear stability 
[73j . Since S'W provides a norm preserved by the lin- 
earized equations, formal stability implies linear stability 
(which in turn implies spectral stability). However, for in- 
finite dimensional systems, formal nonlinear stability need 
not imply nonlinear stability. 

The critical points of energy at fixed mass are given 
by SyV — pSM ~ where is a Lagrange multiplier. This 
leads to u = and to the condition of hydrostatic equi- 
librium 



Vp + pV<P = 0. 



(163) 



Therefore, the critical points of energy at fixed mass are 
steady states of the Euler equation (|150p - ()152|) . For a 
barotropic gas p = Pip), the condition of hydrostatic 
equilibrium (|163p leads to p = p{'P) and p = p{'P) with 
p'{(P) = -p{^)- This leads to the identity 



P'{P) 



p'm' 



(164) 



Since p'{p) = > 0, this relation implies that p'(^) < 0. 

A critical point of energy at fixed mass is a (local) 
minimum iff 



5^W 



J ^ ^^^^^ '^'^^IJ ^^^^ '^''^^ ^^^^^ 



In astrophysics, this is called the Chandrasekhar energy 
principle |21] . 
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for all perturbations 6p that conserve mass. Using Eq. 
(|164|) . this can be rewritten 



dr 



5pS<Pdr>0 (166) 



for all perturbations Sp that conserve mass. We are led 
therefore to considering the eigenvalue equation 



p'{p)Sp + p5<l> = X5p, 



(167) 



with 6'P = u* 6p, for perturbations such that J Spdr = 0. 
For these eigenmodes, we have 



-A 



dr. 



(168) 



Therefore, a critical point of energy at fixed mass is a 
(local) minimum iff all the eigenvalues of the equation 
(|167p are positive. 

Let us specifically consider spatially homogeneous sys- 
tems. Solving Eq. ()167p in Fourier space and using = 
{2'KYu{k)5p since the integral in Eq. (|152p is a convolu- 
tion, we find that the eigenvalues are given by 



A(fc) = + {2TiYu{k)p. 



(169) 



Therefore, a spatially homogeneous system is a (local) 
minimum of energy at fixed mass iff 



cl + {2Trf u{k)p > 0, Vfc. 



(170) 



When condition ()170p is fulfilled, we conclude that the 
system is formally nonlinearly stable with respect to the 
barotropic Euler equations (|150|) - (|152p . We see that, for 
spatially homogeneous systems, the condition of formal 
stability (|170p coincides with the condition of linear sta- 
bility (|154p . On the other hand, when condition (|170p is 
not fulfilled the system is linearly unstable according to 
Sec. 18.11 In conclusion, for spatially homogeneous systems 
described by the barotropic Euler equation, formal nonlin- 
ear stability coincides with linear stability[ 



8.3 Linear stability for the Vlasov equation 

We consider the dynamics of a kinetic system described 
by the Vlasov equation 



df df df ^ 

— + V- — -V^- — = 0, 
dt dr dv 



<?(r,t)= / u{\r~r'\)p{r',t)dv', 



(171) 



(172) 



where u(|r — r'|) is an arbitrary binary potential of in- 
teraction. Clearly, a distribution function / = /(v) that 



This result may remain true for inhomogeneous systems. In 
particular, it is shown in [88] that the neutral modes of linear 
and formal stability for barotropic stars are determined by the 
same condition. 



depends only on the velocity is a stationary solution of the 
Vlasov equation with = pU. We want to investigate the 
linear dynamical stability of this spatially homogeneous 
distribution. Linearizing the Vlasov equation around this 
steady state and using Laplace- Fourier transform (writing 
the perturbations in the form i5/kaj ~ S^Wu> ~ e*''^ ''^'^*^), 
we obtain the dispersion relation : 



, k- ^ 

e{k,uj) = 1 - (27r)''?i(k) / ^ 



k • V — I 



-dv = 0. (173) 



If we take the wavevector k along the z-axis 0, the dis- 
persion relation becomes 



e(fc,w) = 1 - {2Trf u{k) 



kVz — Li 



C f^'Jz 



■dvxdvydvz — 0, 



(174) 



where the integral must be performed along the Landau 
contour C. We define g{vz) = / fdvxdvy. In the following, 
we shall note v instead of and / instead of g. With these 
conventions, the dispersion relation p74p can be rewritten 



e{k,uj) = I- {2Tr fu{k) f I^dv = 0. 



(175) 



We are led therefore to a one-dimensional problem like 
the one investigated in the first part of this paper for the 
HMF model. 

For uJi = 0, the real and imaginary parts of the di- 
electric function e{k,ujr) ~ er{k,uJr) + iei{k,LUr) are given 

by 

er{k,ujr) = 1 - {2TTyu{k)P / ^ ^ dv , (176) 

■J -oo ^ IT 



e,{k,u;r) = -n{2nfu{k)r 



(177) 



The condition of marginal stability corresponds to 
£r(k,ujr) ~ ei{k,LUr) = 0. Thc condition ei{k,uJr) = 
is equivalent to /' (^) = so that ujr/k = vq is equal 
to a velocity where f{v) is extremum (/'(wo) = 0). The 
condition er{k,ujr = wo) = then gives 



l-{2nfu{k) / i-illdi; = 0. 

J-oo V-Vo 



(178) 



This determines the wavenumber(s) k^ corresponding to 
marginal stability. 

On the other hand, as a direct consequence of the 
Nyquist criterion, if /(w) has a unique maximum at v ^ vq, 
then the DF is linearly stable with respect to a perturba- 
tion with wavenumber k iff er{k,u!r ~ vq) > i.e. iff 



1 - {27rf u{k) 



..... /"°°m,„>o. 

oo V-Vq 



(179) 



If the distribution function is isotropic, there is no restric- 
tion in making this choice. 
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This stability criterion generalizes Eq. ([27|) to the case of 
an arbitrary potential of interaction. If u{k) > (repulsive 
potential), a single- humped distribution is always stable. 
This is the case in particular for a Coulombian plasma 
[72] . If u{k) < (attractive potential), the distribution is 
stable to all wavenumbers if 



dv 



> {2Tr)^p\uik)\r 



(180) 



and unstable otherwise. In that last case, the range of 
unstable wavenumbers is determined by the converse of 
inequality (|179[) . For the gravitational interaction in d = 3, 
the stability criterion (|179|) becomes 



1 + 



AttG 

"fc2~ 



+ 00 



f'iv) 
V -vo 



dv > 0. 



(181) 



The system is always unstable for sufficiently small 
wavenumbers: 



k <kj ^ 



1/2 



(182) 



where kj is the Jeans wavenumber for a stellar system. We 
shall make the connexion between the stability of a kinetic 
system and the stability of the corresponding barotropic 
gas in the next section. In particular, using identity (|204p . 
we will show that Eqs. (|179p and (jlSOp are equivalent to 
Eqs. (HMD and (fTTTl) . 



8.4 Formal stability for the Vlasov equation 

The Vlasov equation conserves the mass A/ = J / drdv, 
the energy E = J dvdv + 5 / p^dr and the Casimirs 
Ih = J h(f) drdv where h is an arbitrary continuous func- 
tion. Let us introduce the functionals S = — J C{f ) drdv 
where C is an arbitrary convex function (C" > 0). They 
will be called "pseudo entropies". We also introduce the 
"pseudo free energies" J ^ S - fiE oy F ^ E - TS imd 
the "pseudo grand potential" G = S — pE — aM where 
(3 = 1/T and a are some constants. Since E, M, Ih, S, 
J, F and G are conserved by the Vlasov equation, the 
following optimization principles Ell 



min {E[f] 



all the Casimirs Ih}, 



(183) 



mm {£;[/] I M[f]^M, S[f] = S}, (184) 



max{5[/] I M[f] = AI, E[f]=E}, (185) 



max{J[/] I M[f]=M}, 



(186) 



^"^ A review of these variational principles will be given in a 
forthcoming paper. 



mm {F[f] 



max {G[/]}, 



(187) 



(188) 



determine steady states of the Vlasov equation that are 
nonlinearly dynamically stable. The first criterion is the 
most refined stability criterion that has been introduced 
in the literature [3- The other criteria provide only suffi- 
cient conditions of dynamical stability. Since the solution 
of an optimization problem is always the solution of a more 
constrained optimization problem, we have the relations 



(fT88)i ^ (HHll) ^ (HMl) ^ (HMl) ^ (mi ^ (fT83l) . 



(189) 



We note that ()185p is similar to a condition of microcanon- 
ical (two constraints) stability, p86p or (jl87p is similar 
to a condition of canonical (one constraint) stability, and 
()188p is similar to a condition of grand canonical (no con- 
straint) stability in thermodynamics 15 88]. Therefore, 
the inclusion (fT88)) =^> (fT87)) 44> (fT86| ^ (fTsfe)) is similar 
to the fact that "grand canonical stability implies canon- 
ical stability which implies microcanonical stability" in 
thermodynamics. However, the reciprocal is wrong in case 
of "ensemble inequivalence" that is generic for systems 
with long-range interactions |12I17I18] . In the following, 
we shall mainly be interested by formal dynamical sta- 
bility [73] . We shall consider small perturbations and de- 
termine under which conditions the first variations of the 
functional vanish (critical point) and the second variations 
are positive definite (for problems p83]) . ([T84]l and ([T87]) ) 
or negative definite (for problems ([T85]) . pl6]) and pl8]) ). 

Critical points: considering problem (|183p . we find that 
a DF is a critical point of energy for symplectic perturba- 
tions (i.e perturbations that conserve all the Casimirs) iff 
/(r, v) is a steady state of the Vlasov equation [76|77j . On 
the other hand, the critical points of problems (|184p - (|188p 
lead to DF of the form / = /(e) with /'(e) < 0, i.e. the 
DF depends only on the individual energy e = v'^ /2 + <l>{r) 
and is monotonically decreasing [88] . 

Formal stability: considering problem (|183p . a critical 
point of energy for symplectic perturbations is a (local) 
minimum iff 



/'(^) 



drdv 



6 f 5$ drdv > 0, 



(190) 



for all perturbations Sf that conserve the energy and all 
the Casimirs at first order SE = SIh = 0. Considering 
problems (|184p and (|185p . a critical point of energy (resp. 



In astrophysics, this is called the Antonov energy principle 
[21j . Recently, developing these arguments, Lemon et al. [78j 
have proven for the first time that all DF / = /(e) with /'(e) < 
are nonlinearly stable with respect to the Vlasov-Poisson 
system. 

We could also introduce a grand microcanonical ensemble 
corresponding to the maximization oi K = S — aM at fixed 
energy [15] . 
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pscudo entropy) at fixed mass and pseudo entropy (resp. 
energy) is a (local) minimum (resp. maximum) iff inequal- 
ity (|190[) is satisfied all perturbations Sf that conserve 
mass and energy at first order SM = SE = 0. Consider- 
ing problems (|186|) and p87p . a critical point of pscudo 
free energy J (resp. pseudo free energy F) at fixed mass 
is a (local) maximum (resp. minimum) if inequality (|190[1 
is satisfied for all perturbations Sf that conserve mass 
SM ~ 0. Considering problem (|188p . a critical point of 
pseudo grand potential G is a (local) maximum if inequal- 
ity (|190p is satisfied for all perturbations. Problem (|183p 
is the most refined stability criterion because it tells that, 
in order to settle the dynamical stability of the system, 
we just need considering symplectic (i.e. dynamically ac- 
cessible) perturbations. Of course, if inequality p90p is 
satisfied by a larger class of perturbations, as implied 
by problems (|184p - (|188p . the system will be stable a for- 
tiori. Therefore, we have the implications (|189p . Problems 
(|184p - (|188p provide sufhcient (but not necessary) condi- 
tions of dynamical stability. A steady state can be stable 
according to (|183p while it does not satisfy (|184p - (|185p . 
or it can be stable according to (|184p - (|185p while it does 
not satisfy (|186p - (|187p . or it can be stable according to 
(fT86l) - (fT87l) while it does not satisfy (fTM]) . Therefore, the 
stability criterion (|183p is more refined than (|184p - (|185p 
which is itself more refined than (|186p - (|187p which is itself 
more refined than (fTSH]) . 

Let us consider the minimization problem 

nun{F[/] I M[/] = A/}, (191) 

in more detail. As we have seen, it provides a sufficient 
condition of dynamical stability for the Vlasov equation. 
The critical points of pscudo free energy F at fixed mass 
M are given by 

SF + aTSM = 0, (192) 
where a is a Lagrange multiplier. This yields 

C\f)^-f3e~a. (193) 

Since C is convex this relation can be reversed to give 

f = F{(3e + a), (194) 

where F{x) = (C")"^(-a;). Differentiating Eq. p93)) . we 
obtain 

C"if) = (195) 

so that /'(e) is monotonic. We shall assume that /'(e) < 
so that the pseudo temperature T is positive. A critical 
point of pseudo free energy at fixed mass is a (local) min- 
imum iff 

S^F='^ J C"{f){Sffdrdv+^ J Sp5^dr>0 (196) 

for all perturbations that conserve mass. Using Eq. (|195p . 
this can be rewritten 

S'-F^-^ I ^^drd^ + ^ I SpS^dr>0 (197) 



for all perturbations that conserve mass. 

For any distribution / = /(r, v), the density and the 
pressure are given hy p ~ J f dv and P = ^ / fv^ dv. If 
/ = /(e), we have p = p{<P) and p = p{'P). Eliminating 
^(r) between these expressions we find that the pressure 
is a function of the density: p = p{p). Therefore, to any 
kinetic system described by a distribution function / = 
/(e); we can associate a barotropic gas with an equation 
of state p(p) that is entirely specified by the hmction C(/). 
We now introduce the functional 

F[p] = j Pf^-JT- dp'dr +\ j p1> dr, (198) 

and consider the minimization problem 

min {F[p] I M[p] = M}. (199) 
p 

It can be shown |88I69| that 

(|T9T]l ^ ((1991). (200) 

Therefore, / is a minimum of F[f] at fixed mass iff p is a 
minimum of F[p] at fixed mass. As a result, (|199p provides 
a sufficient condition of nonlinear dynamical stability for 
the Vlasov equation. However, it is not a necessary condi- 
tion of nonlinear dynamical stability. For example, a DF 
can be nonlinearly stable according to (|183p or (|185p al- 
though it does not satisfy (|19ip or p99p . This is similar to 
"ensemble inequivalence" in thermodynamics [12ll7ll8j . 
It is also obvious that 

(fWll ^ (Hm). (201) 

Therefore, if p is nonlinearly stable with respect to the 
Euler equation then it solves (|162p and (|199p . According 
to (|200p . the corresponding distribution function / solves 
(|19ip so it is nonlinearly stable with respect to the Vlasov 
equation. Therefore, a distribution function f = /(e) with 
/'(e) < is nonlinearly stable with respect to the Vlasov 
equation if the corresponding density profile p is nonlin- 
early stable with respect to the Euler equation. In astro- 
physics, this corresponds to the (nonlinear) Antonov first 
law: "a stellar system / = /(e) with /'(e) < is (non- 
linearly) stable with respect to the Vlasov equation if the 
corresponding barotropic star with an equation of state 
p = p{p) is (nonlinearly) stable with respect to the Euler 
equations". However, the reciprocal is wrong in general. 
The distribution function / can be nonlinearly stable ac- 
cording to (|183p or (|185p although it does not satisfy (|19ip 
so that the density profile p does not satisfy (|199p or (|162p . 
We see therefore that the nonlinear Antonov first law is 
similar to the notion of "ensemble inequivalence" in ther- 
modynamics |88j . 

According to the previous discussion, / is a minimum 
of F[f] at fixed mass iff p is a minimum of F[p\ or W[p, u] 
at fixed mass. We can therefore use the results of Sec. 
18.21 In particular, for a spatially homogeneous system, the 
distribution /(v) is a (local) minimum of F[f] at fixed 
mass iff 

cl + {2T:fu{k)p > 0, Vfc, (202) 
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where cj. ~ p'{p) is the velocity of sound in the correspond- 
ing barotropic gas. Now, using Eq. (|164|1 and the fact that 
p = J /(e) dv, we have for any distribution / — /(e): 



p(r) 



P(r) 



p{r) 



J Vz J—OC Vz 



(203) 



where g{vz) = J /(v) dv^dvy. For a spatially homogeneous 
system, we obtain 

(204) 



+ 00 f'(v) 



dv 



where we have noted v for Vz and / for g like in Sec. 
18.31 Substituting this relation in Eq. (|202p . we find that 
a homogeneous distribution is a (local) minimum of F at 
fixed mass iff 



1 - i2Trfu{k) 



+ 00 



/'(«) 



dv > 0, Vfc. 



(205) 



because it is valid for an arbitrary potential of interac- 
tion. Our method exploits the link between the variational 
problems acting on the distribution function / and those 
acting on the spatial density p. We have used the fact 
that a distribution function / is stable with respect to the 
Vlasov equation if the corresponding density profile p is 
stable with respect to the Euler equation. This is similar 
to the Antonov first law in astrophysics and this is related 
to a notion of ensemble inequivalence in thermodynam- 
ics (for spatially homogenous systems, we have ensemble 
equivalence). As a by-product, our approach gives a phys- 
ical interpretation of the stability criterion ()205|1 in terms 
on a condition on the velocity of sound (|202p in the cor- 
responding barotropic gas. Furthermore, our method can 
be extended to spatially inhomogeneous systems (see Ap- 
pendix[F|. Indeed, if all the eigenvalues A of the eigenvalue 
equation (|167p are positive, then p is formally stable with 
respect to the Euler equation, implying by ()20f |) and (j2GG|) 
that the corresponding distribution / is formally stable 
with respect to the Vlasov equation. 



When condition (|205p is fulfilled, we conclude that the 
system is formally nonlinearly stable with respect to the 
Vlasov equation [according to criterion (|187[) ]. We see 
that, for spatially homogeneous systems, the condition of 
formal stabilit y (I205D coincides with the condition of linear 
stability (fT79)) Fl. When condition (|205)) is not fulfilled the 
system is linearly unstable according to Sec. 18.31 There- 



fore, for homogeneous systems, criterion (|187p is a suffi- 
cient and necessary condition of dynamical stability (in 
that case, we have ensemble equivalence). In conclusion, 
for spatially homogeneous systems described by the Vlasov 
equations, formal nonlinear stability coincides with linear 
stability. On the other hand, Eq. (|2G4p shows that the 
integral appearing in Eqs. (|205p and (|179p is directly re- 
lated to the sound velocity in the corresponding barotropic 
gas Ell- This implies that the condition of (linear and for- 
mal) dynamical stability for the Vlasov equation coincide 
with the condition of (linear and formal) dynamical sta- 
bility for the Euler equation. In conclusion, a spatially 
homogeneous system is dynamically stable with respect to 
the Vlasov equation iff the corresponding barotropic gas 
is dynamically stable with respect to the Euler equations. 
This is the proper formulation of the Antonov first law 
for homogeneous systems. In that case, we have "ensem- 
ble equivalence" . This is not true in general for spatially 
inhomogeneous systems. 

Remark: we have obtained the condition of formal sta- 
bility (|205p by a method different from the one devel- 
oped by Yamaguchi et al. [32j . Our result is more general 



2° Since / = /(e) with /'(e) < where e = v'^ /2 + #(r), we 
deduce that, for a spatially homogeneous system {$ = est.), 
f{v) is a symmetric function with a single maximum at vo = 0. 
Therefore, Eq. ([T79|l reduces to Eq. p05|l . 

'^^ This is true for any distribution function, not only for the 
Maxwell and Tsallis distributions. However, the relation (|204|l 
can be checked explicitly for these simple distributions |29j . 



9 Conclusion 



In this paper, we have carried out an exhaustive study of 
the dynamical stability of systems with long-range inter- 
actions described by the Vlasov equation. For illustration, 
we have considered the HMF model and we have treated 
the case of single-humped (Maxwell, Tsallis,...) as well 
as double-humped (symmetric and asymmetric) distribu- 
tions. The stability of the system has been investigated by 
using the Nyquist method. Interestingly, the stability dia- 
grams of these systems are very rich, exhibiting re-entrant 
phases in the case of double- humped distributions. For the 
HMF model, the potential of interaction is truncated to 
one Fourier mode so that only the modes n = ±1 can 
propagate. As a result, the stability of the system is ex- 
pressed as a condition on the temperature T while for 
other systems (plasmas, gravitational systems,...) it is ex- 
pressed as a condition on the wavenumber k (for a given 
value of the temperature T). On the other hand, for the 
HMF model, we can consider attractive as well as repul- 
sive interactions and there exists spatially homogeneous 
distributions in each case allowing the application of the 
Nyquist method. Therefore, by changing the sign of the 
interaction, this simple toy model exhibits features sim- 
ilar to Coulombian or gravitational plasmas. Noticably, 
this is the first time that the Nyquist method is applied 
to a system with attractive interactions without making 
any (Jeans) swindle. Finally, we have shown that the re- 
sults and methods developed for the HMF model are in 
fact applicable to more general potentials of interaction 
giving further interest to our study. This opens the way 
to many numerical and/or analytical investigations that 
will be considered in future works. 
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A The function Wr{z) 

Let us study the function 

Wr{z) ^ 1 - ze 2 / e 2 dx, 
Jo 

for z real. We clearly have 

Wri-Z) ^ Wr{z), 



(206) 



(207) 



so we can restrict ourselves to z > 0. This function is plot- 
ted in Fig. [ST] It starts from 1 for 2 = 0, crosses the a;-axis 
for z = Zc = 1.307, achieves a minimum value —0.2847 for 
z ~ 2.124 and tends to for z — > +oo. In order to deter- 
mine its asymptotic behavior, let us consider the integral 



1=1 e dx. 

'0 



(208) 



It can be rewritten 



,2_^2 

dx. 



(209) 



With the change of variables y — \/z^ ~ x'^, we get 

y 



/ = -e 2 

z 



For z +00, we obtain 



1 II 



:dy. 



(210) 



r 1 r , f. ^y^ 

I=-e^ I dye-^y(l + f-^ + ^ 



so that finally 



1 .2 / 1 3 

I=-e- 1 + — + — 

z \ z^ z* 

Therefore, we get 

1 3 

z 

On the other hand, for z — > 0, we have 

Wriz) = 1 - - .. 



(211) 



(212) 



Wriz) 



(z^+oo). (213) 



(z ^ 0). (214) 



B The Tsallis and Cauchy distributions 

The Tsallis distributions are usually written in the form 
[83]: 



l-('?-l) 



26 



(215) 



where Cq is a normalization constant (see below). The 
parameter 9 > has the dimension of a temperature but 
it is not the kinetic temperature T = (w^). To arrive at 




Fig. 61. The function Wr{z) for z real. 



the form (|47p . we can either follow the steps described in 
|29|83j or directly proceed as follows. First, we introduce 
notations similar to those used in astrophysics |21j and 
define the polytropic index n through the relation l/(g — 
1) = n — 1/2. Equation (|215p can then be rewritten 



/ = A 



1 - 



1 



129 



n-l/2 



(216) 



This distribution is normalizablc provided that n > 1/2 
or n < 0. Using p — J f dv, we get 



r{n + i) 



r(7i H- 1/2) (n- 1/2)1/2' 

r(l/2- n) 
r(-n)[-(n- l/2)]i/2' 



n>-, (217) 



n < 0. 



(218) 



According to criterion (|28p . the Tsallis distributions (|216p 
are linearly stable with respect to the Vlasov equation iff 



> 



kM 



71-1/2 Att 



(219) 



On the other hand, the kinetic temperature (variance) 
T = (u^) exists if n > 1/2 or n < —1. In that case, we 
have 



T = 



(220) 



Substituting this relation in Eq. (|216p yields expression 
(|47p of the distribution function where T has a clear phys- 
ical meaning as being the kinetic temperature (or the ve- 
locity dispersion). 

The Tsallis distributions with index — 1 < ti < have 
not been considered in Sec.[3]because they have infinite ki- 
netic energy K = N{v'^)/2, so they are not physical. Yet, 
they can be studied at a mathematical level. For these dis- 
tributions, the stability criterion is given by Eq. (|219p and 
the Nyquist curves are similar to those reported in Sec. [3] 
for n < — 1. Of particular interest is the Cauchy distribu- 
tion corresponding to n = —1/2. The Cauchy distribution 
is 



TT V29 1 
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and the corresponding dieleetric function ^ can be writ- 
ten 



kM 



An^e Jc (l+x2)2(a; 



/2e 



dx. (222) 



The integral can be computed analytically by complet- 
ing the domain of integration by a large semi-circle in 
the lower-half plane of the complex plane and using the 
Cauchy residue theorem for a function with a pole of order 
2 at X = — i. This yields 



1 



kM 



1 



The condition e{Lo) ~ gives the complex pulsations 



(223) 



(224) 



For the attractive HMF model (see Sec. [2]), the pulsations 
are purely imaginary. The mode lo+ is always damped 
(stable). The mode w_ is damped (stable) if 6* > while 
it grows exponentially with time (unstable) ii 6 < 9c with 



kM 
"87' 



(225) 



in agreement with Eq. ()219p . On the other hand, if we 
take LOj = 0, we find that 



kM 1 



26 



ei{ujr) = - 



8716* 



kM 



1 + ^ 



V2e 



iTT0 



(226) 



(227) 



The corresponding Nyquist curve is represented in Fig. 
[62l It is very similar to the one obtained for other values 
of n < but, for the Cauchy distribution, it is obtained 
analytically. 

The case of the repulsive HMF model (see Sec. [7]) is 
obtained by changing +k in —k in the preceding relations. 
The condition e(uj) ~ gives the complex pulsations 



kM 
47r 



iV20 



These modes are always stable. They evolve with the 
(proper) pulsation ujp = ± ^/kM/^A-K) (see Appendix |E| 
independent on the temperature and they are damped at 
a rate 7 = ^/W. For cji = 0, the real and imaginary parts 
of the dielectric function are given by Eqs. (|226p and (|227p 
with k replaced hy —k. The corresponding Nyquist curve 
is represented in Fig. [62l 
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Fig. 62. Nyquist curve (analytical) for the Cauchy distribution 
for the attractive and repulsive HMF models. 



C The function W^''\z) 

Let us consider the function 



/27r n 



(229) 



for z real and n < — 1 or n > 1/2. We clearly have 

^^(")(-z) Wf^{z). For n < -1 and n > 5/2, this 
function is defined for any z and its derivative is continu- 
ous. For 3/2 < n < 5/2, this function is defined for any z 
but its derivative is discontinuous at 2: = y/2{n + 1). For 
1/2 < n < 3/2, this function diverges at z — \/2{n + 1). 
This leads to the behaviors represented in Figs. [63] and 
M For n > 3/2, we define (n = -W^''\^/2{n + I)) or, 
explicitly. 



/2tt n 



n 

2 



2(n+l)J 



ln-3/2 



y/2{n+l) V2(n+T) 



-dx. 



(230) 



Note that the principal value is not needed in that case. 
Note also that Cn ^ 0. This inequality can be obtained 
from Eq. ([25]) by taking Vext = \/2{n + 1) and recalling 



that f{vext) = 0. Setting y = x/y^2{n+l), Eq. (^50)1 can 
be rewritten 



C.~^U4)v/2(;rTT)/„, (231) 



(228) with 



+1 



y 



{l^y)^-5/^l+y)n-3/2^y^ (232) 



It turns out that this integral can be computed analyti- 
cally yielding 



/„ = ^r{n - 3/2) 



4"r(n-hl/2) 40F 



r(2n-l) r(n-l) 



(233) 
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Fig. 63. The function Wr'{z) for n < -1 and n > 3/2 




Fig. 64. The function W^"'{z) for 1/2 < n < 3/2 
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On the other hand, for z E [— VS, VS], we have 



(3/2) 



2^/5 ' 



(235) 



while W.l (z) = otherwise. Using these results, the 
Nyquist curve in Fig. [5] can be obtained analytically. 

Finally, for the index n — 3/2 and for z = ij where 
7 is a real number, the function (|54p can be calculated 
analytically. For 7 > 0, using Eq. (fTO|) . we get 



(236) 



W^(3/2)(,^)^l_^arctanf^ 



For 7 < 0, using Eq. (fT^ . we get 



v5 



arctan 



(237) 



Finally, for 7 0, using Eq. (HH), we have W'-^/'^^O) = 0. 
Therefore, looking for a solution of the equation e{n) = 
1 - {v/Vc)W'^^^^{^n) = of the form fi = ifi,, we find 
that for 77 > = 5/3 (unstable), the perturbation grows 
with a rate f2i > given by 



1-^ 

Vc 



V5 



arctan 



' V5 \ 



0, (238) 



while for 77 < ?7c = 5/3 (stable), the perturbation is 
damped with a rate f2i < given by 



1-^ 



I _ ^ ' i arctan [ ^ | + tt 

V5 \ i^v^^^.y ^ 



0. (239) 



Fig. 65. The parameter = -W^r*"'( 1/2(71 + 1)) as a function 
of the polytropic index n > 3/2. 



For n 3/2+, we have the equivalent C„ ^ l/[2(n— 3/2)]. 
For 77, +00, (n 0. This parameter is plotted as a 
function of n in Fig. 1651 

Finally, for the index n = 3/2, the functions ([57)) - (|58p 
can be calculated analytically. For any real z / ±i/5, we 
find 



Finally, for rj ^ rjc = 5/3, we recover the case of marginal 
stability i?; = 0. For 77 +00, we recover Eq. ()276p and 
for 7/ 0, we find 5\/5/(37r773/2). 



D Analogies and differences between the 
HMF model, Coulombian plasmas and 
self-gravitating systems 

For Coulombian plasmas, the potential of interaction is de- 
termined by the Poisson equation Au = — 47r(e^/777^)(5(x) 
yielding 



(2Trf u{k) 



2^/5 



In 



V5- 



V5- 



(234) 



We introduce the plasma frequency 

= 2~' 

and the Debye wavenumber 



> 0. 



(240) 



(241) 



(242) 
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The dispersion relation (|175p can be written 
(I ^ 1 47re2 [ f{v) 



(243) 



Comparing this expression with Eq. (|127[) . we make the 
link between Coulombian plasmas and the repulsive HMF 
model by setting 



= -, /3m = /?, -=u, -^=^p, (244) 



where the tilde variables refer to the HMF model and the 
un-tilde variables to the plasmas. Then, we find that the 
dimensionless parameters for the plasmas are 



fr2 

fc2 ■ 



(245) 



The results obtained in Sec.[7]for the repulsive HMF model 
can be applied to Coulombian plasmas provided that we 
interpret as a normalized squared wavelength and as a 
normalized pulsation according to Eq. (|245p . In particular, 
single-humped distributions are stable to all wavelengths 
[72] and the Nyquist curves for Maxwell and Tsallis dis- 
tributions are given in Sec. 17.1.11 By contrast, the phase 
diagrams of the double- humped distribution (in Sec. 17. 1.2)) 
are not directly applicable to plasmas because, according 
to Eq. (|245p . the dimensionless parameter a cx l/y^ de- 
pends on k. Therefore, this is not the correct dimensionless 
separation. For plasmas, the correct dimensionless separa- 
tion (which must be independent on the wavenumber k) is 
y = T]a'^ ~ Pmv'^ instead of a. Therefore, the correct sta- 
bility diagram corresponds to the curve(s) ric{y) obtained 
from Eqs. ()137|) - (|138p . This will be discussed in more de- 
tail in a specific paper. 

For self-gravitating systems, the potential of interac- 
tion is determined by Au — 47rG'(5(x) yielding 



{2TTfu{k) 



AttG 

1^ 



< 0. 



We introduce the Jeans wavenumber 

kj = ATrGPpm. 
The dispersion relation (jl75p can be written 



e(fc,w) 



/'(«) 



dv. 



(246) 



(247) 



(248) 



Comparing this expression with Eq. ^ , we make the link 
with the attractive HMF model by setting 



47rG 



l3m = (3, 



k 



(249) 



where the tilde variables refer the HMF model and the un- 
tilde variables to self-gravitating systems. Then, we find 
that the dimensionless parameters for self-gravitating sys- 
tems are 



77 



[2 



(250) 



The results obtained in Sec. [5] for the attractive HMF 
model can be applied to (uniform) self-gravitating sys- 
tems provided that we interpret 77 as a normalized squared 
wavelength and 12 as a normalized pulsation according to 
Eq. (|250p . In particular, the Maxwell distribution is stable 
for k > kj (i.e. 77 < 1) and unstable for k < kj (i.e. 77 > 1) 
and the Tsallis distributions are stable for k > (—TrY^^k j 

V n+l ^ 

(i.e. 77 < 77c) and unstable for k < {i^^Y^^kj (i.e. 77 > 77c). 
The Nyquist curves for the Maxwell and Tsallis distribu- 
tions are given in Sees. [3] and [D By contrast, the phase 
diagrams of the double- humped distribution (in Sees. 12.61 
and l2.7p arc not directly applicable to self-gravitating sys- 
tems because, according to Eq. (|250p . the dimensionless 
parameter a (x l/^/rj depends on k. Therefore, this is not 
the correct dimensionless separation. For self-gravitating 
systems, the correct dimensionless separation (which must 
be independent on the wavenumber fc) is y = 770^ = (imv^ 
instead of a. Therefore, the correct stability diagram cor- 
responds to the curve(s) rjciy) obtained from Eqs. (|llip - 
()112p . This will be discussed in more detail in a specific 
paper. 



E Landau damping and Landau growth for 
the repulsive and attractive HMF models 

Like for Coulombian plasmas, the repulsive HMF model 
leads to Landau damping at low temperatures. On the 
other hand, the attractive HMF model can lead to "Lan- 
dau growth" associated to the phase transition (collapse) 
at low temperatures. 



E.l The limit T ^ 

Let us first consider the limit T ^ Q for the repulsive 
HMF model. The case T = corresponds to f = p S{v). 
In that case, we obtain e(a>) = 1 — The condition 
e(w) = yields uj = ±(fcp/2)i/2 
proper pulsation 



±ujp where is the 



kp 
T 



kM 
An ' 



(251) 



This is the equivalent of the plasma pulsation in plasma 
physics. Therefore, for T = the system is stable and 
the perturbation oscillates with the pulsation LUp without 
attenuation. To determine the damping rate 7 — —LUi for 
T <C 1, we shall solve the equation (|127p for the complex 
pulsation u = ujr+iuJi by assuming uji ^ lUr (we adapt the 
method of plasma physics given in [72| ) . Expanding the 
dielectric function in Taylor series as e{uj) = e{ujr + ~ 
e{LUr) + ie'{ujr)uji, the condition eiuj) = gives to leading 
order 



r{uJr) = 0, 



i{uJr) + e'r{^r)(^i = 0, 



(252) 



(253) 
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where er(wr) and ei{ujr) are defined in Eqs. (|128|1 and 
(|129|) . Integrating by parts, Eq. ()128p can be rewritten 



the form of the distribution. For the MaxwcUian distribu- 
tion we find that 



ericor) = 1 - |P / °° dv. (254) 

For w/o^r ^ 1, we get 



e,K) = 1 - -il ( 1 + ^ + 



2tj, 



(255) 



If the distribution is symmetrical with respect to the origin 
and if the variance is finite, we obtain 



kp f 3(w2\ 



and 



(256) 



(257) 



Substituting these relations in Eqs. (|252p - (|253p . we find 
that 



i-^ri+^i=o. 



fc , kp 

2 Ci),, 



(258) 



(259) 



This expansion is valid for T = (u^) ^ 0. For T = 0, Eq. 
(|258p gives Wr = Wp. Then, to next order, Eq. (|258p gives 



uj-;. = + 3T - 



(T^O). 



(260) 



This asymptotic expansion for T — *■ also gives the ex- 
act pulsation (valid at any temperature) for the water-bag 
distribution (see Sec. 17.1.1]) . Using the dimensionless vari- 
ables introduced in Sec. [31 Eq. (|260p can be rewritten W\ 



{r\ +oo). 



(261) 



On the other hand, to leading order when T (or 
r\ +oo), Eq. ([259]) gives the damping rate 



p fi 



2p 



(262) 



When f'{u!p) < 0, we obtain the expression of the Landau 
damping for the repulsive HMF model. We stress that the 
relation (|260p is the same for any distribution. By con- 
trast, the damping coefficient (|262|) strongly depends on 



■^■^ Let us recall that this solution only exists for the repulsive 
HMF model. For the attractive HMF model, there is collapse 
when T —> (see below) . 

A more accurate expression is obtained by replacing f'{cjp) 

by f'{iOr). 



^3/2g-r,/2^ 



(263) 



For the asymmetric double- humped distribution (|93p . we 
get 



1 /^^3/2 



(1 - a)e-^(i-'^)' 



-HZ\(1 + a)e-*(i+")' 



For the Tsallis distributions ([47]) . we obtain 



2 / n + 1 



(264) 



-| ji-3/2 



1 - 



2(n + l) 



(265) 



For n ±oo, we recover Eq. (|263[) . For n > 1/2, the 
damping is zero if LOp > i.e. rj > 2{n + 1). For fixed n 
in the range ] — oo, — 1[ and i] +oo, Eq. (|265p becomes 



1 

n — — 



2) [-(n + l)r 



-tT, (266) 



so the damping rate behaves like fli cx — I/t^I"' for 77 
+00. Let us consider for illustration the index n = —3/2 
corresponding to the distribution 



/ 



2p 



ttVT (1 + ^)' 



In that case, we have 



l3/2- 



(267) 



(268) 



Using the correspondance between Coulombian plas- 
mas and the repulsive HMF model (see Appendix [P]) . we 
find that the equivalent of Eqs. (|260p and (|262p are 



=u? + 3Tk^ 



2pk'' 



(^) . (269) 



These expressions are valid for k/kn 0. For the 
Maxwellian distribution, we get 



/ TT ( k£i 



(270) 



For the Tsallis distributions, we get 



ko 
k 



1^ 1 

^ \ n — - 



2 1 n+1 



1 - 



2{n+l)k^ 



-3/2 



(271) 
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For n — > ztoo, we recover Eq. (|270p . For n > 1/2, the 
damping is zero if ojp/fc > Vm, i-e- iik/kj < 1/ y^2(n + 1). 
For fixed n in the interval ] — oo, — 1[ and k/k]j — > 0, we 
obtain 



ITT 1 



2"- 2 



1 



2 y f- 



1 



(n + l)]' 



kp 



2n 



(272) 



so the damping behaves like uji oc for k/ko 0. 

For the distribution (|267p corresponding to n = —3/2, we 
have 



4CJ; 
"P" 



(273) 



Let us now consider the limit T — > for the attractive 
HMF model. The case T = corresponds to / = p (5(f). 
In that case, we obtain e{uj) = 1 + = yielding 
Lo = ±i{kp/2)^/'^. Therefore, for T = the system is 
unstable and the perturbation grows with a growth rate 
oji = {kM/AnY^^. Let us now consider the case T <C 1. 
We anticipate a solution with cOr = and uJi > 0. In that 
case, the perturbation grows (tUi > 0) without oscillating. 
The growth rate is given by Eq. ([30)1 . Integrating Eq. (|30p 
by parts, we obtain 



k 

'--2 



+00 



2\2 



dv = 0. 



(274) 



Expanding the integrand in powers of v/tOi <JC 1 for T 0, 
we find that 



kM 
An 



3T 



(T^O), 



(275) 



with T — {v^). This asymptotic expansion for T — > 
also gives the exact pulsation (valid at any temperature 
T < Tc = kM/\2iT) for the water-bag distribution (see 
Sec. H]). Introducing dimensionless variables, we get 



V 



[rj +00). 



(276) 



For the Maxwellian distribution, we recover Eq. (136) of 



Separating real and imaginary parts, we obtain two tran- 
scendant equations 



Re [ikTrf'{LUr + = 1, 



Im [if'{ujr + i^i)] = 0, 



(279) 



(280) 



which strongly depend on the form of the distribution. 
For the Maxwellian ([34|) . they can be rewritten to leading 
order in the limit uji/cur 3> 1 as 

/ o \ 1/2 

fcTT f ^ j pPujie^"^^ cos(/5wrwO = 1, (281) 



Sm{f3LUrLUi) — 0. 



(282) 



Equation (|282p implies f3oJrOJi = rmr. Eq. (|28ip will have a 
solution provided that m is even. Let us take Pojr^i = tt. 
Then, Eq. (|28T|) gives 



(283) 



which determines coi. By a graphical construction, it is 
easy to see that \uJi\ is an increasing function of T. For 
T +00, we find the asymptotic behavior 



V2rinT, 



T 



2\nT 



(284) 



Since LOi/cUr ^ InT ^ +00, our basic assumption is sat- 
isfied. We also note that < in the repulsive case. 
Introducing dimensionless variables, we have for 77 — > 0: 



f2, 



/2|lnr/| 



(285) 



Using the correspondance (|245p between the repulsive 
HMF model and Coulombian plasmas, we recover the re- 
sults of plasma physics [T^ for k/ko ^ +00: 



2ujf 
ko 



fcv In fc. 



2k D \/hik 



(286) 



E.2 The limit T -> +00 

Let us consider the limit T — > -l-oo for the repulsive HMF 
model (we adapt the method of plasma physics given in 
[7D]1. We look for a solution of the equation e{uj) = of the 
form Lu = LUr + iLUi with LUi < (damping) and LUi/ur ^ 1. 
Using Eq. ([12]) for cj^ < , Eq. (fT27l) can be written 



1 



/'(«) 



dv — ikiT f' {bj) = 0. 



(277) 



If / decreases exponentially rapidly for large w , then for 
Ui/iOr 3> 1, the foregoing equation reduces to 



1 - ikirfiuj) = 0. 



(278) 



Let us consider the limit T —f +00 for the attractive 
HMF model. We again look for a solution of the equation 
e{Lu) = of the form uj = ujr + iloi with tj^ < (damp- 
ing) and LOi/iDr 3> 1. Repeating the procedure followed 
previously, we obtain two transcendant equations 



Re [ikTrf'{uJr + ii^i)] 



-1, 



Im [if'{Ljj. + iuJi)] = 0. 



(287) 



(288) 



For the Maxwellian ([34|l . they can be rewritten to leading 
order in the limit LOi/uOr ^ 1 as 

/ ^ \ 1/2 ^ 
fcTT f ^ j p/3u;,e^"' cos{(3L0ri^i) = -1, (289) 
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sin(/?WrWi) = 0. 



(290) 



Equation (|290|1 implies f3ujr^i ~ mn. Eq. (|289p will have a 
solution provided that m is odd. Let us take uJr = 0. Then, 
Eq. ([289]) reduces to Eq. ([2831) . Therefore, for T +oo, 
we find the asymptotic behavior 



0. 



(291) 



uji = —V2T InT, LUr = 
This returns Eq. (141) of 

F Marginal stability of spatially 
inhomogeneous states 



For the HMF model, the eigenvalue equation associated 
with the minimization problem (|162p can be written (see 
Eq. (77) of EHl) 



de\ p de 



k 

2^ 



e')de' 



-2Xq, 
(292) 



with q{0) = q{2n) = 0, where q{9) = 5p{e')de' is the 
perturbed integrated density This equation is valid for 
homogeneous and inhomogenous equilibrium density pro- 
files p{e). We have S'^W ^ \ j C[q]qde where £ is the 
operator defined by the r.h.s. of Eq. ()292p . Therefore, a 
critical point of energy yy[p, u] at fixed mass is a (local) 
minimum iff all the eigenvalues A are positive. In that 
case, the density profile p{6) is formally nonlinearly stable 
with respect to the Euler equation. On the other hand, 
since (|162p <^ (|199p <^ (|19ip . the corresponding distribu- 
tion function f{9, v) (leading to an equation of state p{p)) 
is a (local) minimum of pseudo free energy F[f] at fixed 
mass, so that it is formally nonlinearly stable with respect 
to the Vlasov equation. Therefore, the positiveness of the 
eigenvalues of Eq. (|292p is a sufficient condition of formal 
nonlinear stability for the Vlasov equation. Unfortunately, 
it does not appear feasible to solve Eq. (|292p analytically 
in the general case. However, it is possible to solve it an- 
alytically for the marginal case A = 0. This criterion will 
give the point at which the series of equilibria of the inho- 
mogeneous equilibrium states becomes unstable. This is 
the first result concerning the stability of inhomogeneous 
states of the HMF model. The eigenvalue equation (|292p 
can be rewritten 



_^ f p'jp) dq 
de\ p de 



A^ 

— bx cos 9 — by sin 6 = — q, (293) 



with 



bx = — / q{0)cosede = — q'{e)smede, 

2tt Jo 27r 



(294) 



Note that the eigenvalue equation associated with the lin- 
ear stability problem (see Eq. (85) of [29]) has a similar form 
as Eq. (|292|) and that they both coincide at the neutral point 
A = 0. 



- I q{e)sm9de = -^ I q'{e)cos9d9, 
Jo 27r Jq 



(295) 



where we have used an integration by parts to obtain the 
second equalities. For A = 0, Eq. (|293p is readily integrated 
to yield 



dq 
d0 



P 



P'{P) 



{bx sin — by cos 9 + C), 



(296) 



where C is a constant of integration. Another integration 
with the boundary condition g(0) = yields 



'Z(^) 



sine' de' 



p 



cose' de' + c 



■de'. (297) 



" lo P'ip) ' ^ Jo P'ip) 

Then, the condition q{2TT) = determines the constant 



C = b,, 



CyuA^osede 

/c 



p'jp) 

de 

p'(p) 



(298) 



where we have assumed that the density profile p{e) is 

gmmetric with respect to = to simplify some terms 
. Substituting Eq. ([296]) into Eq. ([294]) . we find either 
that 6^ = or that 



i = A 

27r 



P'ip) 



sin^ e de. 



(299) 



Substituting Eq. 
by ^ or that 



into Eq. (|295p we find either that 



k 



27r 



p'ip) 



cos^ede 



( _p_ 

\J0 p'ip) 



cos e de 



27r 



f'^ ^ de 

JQ p'ip) 



(300) 



Therefore, if we restrict ourselves to even perturbations 
{by = and bx ^ 0) the condition of marginal stability is 
given by Eq. (|299p and if we restrict ourselves to odd per- 
turbations {bx = and by ^ 0) the condition of marginal 
stability is given by Eq. (|300p . For general perturbations, 
we must satisfy both Eqs. (|299p and (|300p . Note that using 
Eq. (|203p . which becomes for the HMF model 



p'ip) 



dv. 



(301) 



we can express the criteria of marginal stability (|299p 
and (|300p directly in terms of the distribution function 
f{e,v) = f{e). 



Since p — p{$) at equilibrium with $ = cos d + By sin 9 
(where B^;, By are the two components of the magnetization), 
it is always possible to choose the x-axis such that the distri- 
bution profile is even {By = 0). 
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For homogeneous distributions p = const., we imme- 
diately find that Eqs. ([299]) and ((3001) both reduce to the 
condition of marginal stability 



kM 
1^' 



+00 



dv = 0. (302) 



For inhomogeneous isothermal distributions (see Eq. 
(20) of [29]) for which p = pT, using Eqs. (21) and (26) of 
[29] and the relation li (t) = Iq (t) — jh (t) , we find that Eq. 
(|299p is always satisfied. This implies that an inhomoge- 
neous isothermal distribution is always marginally stable 
with respect to even perturbations {by ~ and bx 7^ 0) 
of the form (|297p . for any temperature T. This can be 
checked explicitly. For an isothermal distribution, these 
perturbations can be written 



-PBcose 



1, 



(303) 



and a direct substitution shows that they are solution of 
Eq. (|293p with A = 0. Finally, for inhomogeneous isother- 
mal distributions, using Eqs. (20), (21) and (26) of [25] , 
we find that Eq. ()300|1 is equivalent to 



2 /i(x) 
X loix) 



(304) 



with X = j3B. The only solution of this equation is a; = 
leading to T = = T^. Indeed, it is shown in [SSj 
that the inhomogeneous isothermal distributions are al- 
ways stable (A > 0) so that the point of marginal stability 
corresponds to the point a; = where the branch of in- 
homogeneous states begins. At the temperature T — Tc, 
the branch of homogeneous states becomes unstable and 
bifurcates to the branch of inhomogeneous states (sec- 
ond order phase transition). Therefore, in that case, the 
marginal stability condition for homogeneous and inhomo- 
geneous states coincides. However, this is not necessarily 
the case for more general distributions like for example the 
Lynden-Bell (or Fermi-Dirac) distribution where a first or- 
der phase transition takes place [57|58j. In that case, the 
new stability criteria (|299p and (|300p can give informa- 
tion about the stability/instability of the inhomogeneous 
states. 
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